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FOREWORD 


Study  to  adapt  Supporting  Design  Synthesis  Programs  for  cost  estimating  use  was 
conducted  by  the  Convair  Aerospace  Division  of  General  Dynamics,  San  Diego,' 
California,  under  USAF  Contract  F33615-72-C-2083.  The  contract  titled  "Weapon 
System  Costing  Methodology  for  Aircraft  Airframes  and  Basic  Structures, " was 
initiated  under  Project  1368,  "Advanced  Structures  for  Military  Aerospace  Vehicles," 
Task  136802,  "Structural  Integration  for  Military  Aerospace  Vehicles."  The  work 
was  administered  under  the  direction  of  the  Air  Force  Flight  Dynamics  Laboratory, 
Structures  Division,  Wright-Patterson  Air  Force  Base,  Ohio,  under  the  direction  of 
Mr.  R.  N.  Mueller  (AFFDL/FBS)  as  Project  Engineer. 

This  report  covers  work  conducted  from  July,  1972  to  November,  1973  and  was 
submitted  by  the  author  in  December,  1973,  under  General  Dynamics  Report  CASD- 
AFS-73-001  as  an  Interim  Technical  Report.  The  principal  author  and  proiect 
leader  on  this  program  is  Mr  R.  E.  Kenyon  of  Convair  Aerospace.  Others  who 
contributed  to  these  studies  include  Messrs.  B.  II.  Oman  and  W.  D.  Honeycutt, 

Mass  Properties;  and  L.  M.  Peterson  and  J.  R,  Hotz,  Structural  Analysis. 


This  report  has  been  reviewed  and  is  approved  for  publication. 
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ABSTRACT 


This  volume  describes  the  supporting  programs  used  in  conjunction  with  a cost  esti- 
mating program  to  provide  a trade  study  cost  estimating  technique  for  aerodynamic 
surfaces.  The  supporting  programs  for  the  purpose  of  this  discussion  are  defined  as 
a structural  synthesis  program  and  a secondary  structure  synthesis  program.  The 
structural  synthesis  program  is  used  for  the  analysis  of  primary  structure  and  is 
called  APAS  (Automated  Program  for  Aerospace -Vehicle  Synthesis).  The  secondary 
structure  synthesis  program  estimates  geometry  and  weights,  and  performs  parts 
definition  for  the  aerodynamic  surface  leading  edge,  trailing  edge,  and  tip  com- 
ponents. The  cost  estimating  adaptation  is  derived  from  programs  originally  deve- 
loped under  Independent  Research  and  Development. 

An  operating  interface  between  cost  estimating  procedures  and  the  synthesis  pro- 
grams has  been  developed.  Dimensional  and  weight  data  from  the  synthesis  pro- 
gram outputs  a e used  as  explanatory  variables  in  a series  of  cost  estimating 
relationships.  The  two  synthesis  programs  also  interface  with  and  require  inputs 
from  pre-design  activities  that  may  also  involve  computerized  programs.  The 
various  capabilities  and  input  requirements  of  the  synthesis  programs  are  described. 
References  are  provided  for  those  desiring  to  pursue  further  the  technical  descrip- 
tions. 
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SECTION  I 
INTRODUCTION 

This  volume  describes  the  supporting  programs  used  in  conjunction  with  the  cost 
estimating  program  to  provide  a trade  study  cost  estimating  technique  for  aero- 
dynamic surfaces.  The  flow  diagram  of  this  estimating  method,  repeated  from 
Volume  1,  is  shown  as  l-igure  1.  The  supporting  programs  for  the  purpose  of  this 
discussion  are  defined  as  the  structural  synthesis  program  and  the  secondary 
structure  synthesis  program.  The  structural  synthesis  program  is  used  for  the 
analysis  of  primary  structure  and  is  called  APAS  (Automated  Program  for  Aerospace 
Vehicle  Synthesis).  The  secondary  structure  synthesis  program  estimates  geometry 
and  weights,  and  performs  parts  definition  for  aerodynamic  surface  leading  edge, 
trailing  edge,  and  tip  components.  The  operating  interface  between  cost  estimating 
procedures  and  the  synthesis  programs  is  illustrated  in  Figure  2.  These  two  pro- 
grams also  interface  with  and  require  inputs  from  pre-design  activities  that  may 
also  involve  computerized  programs.  This  interface  is  shown  in  Figure  3.  Although 
the  trade  study  cost  estimating  method  can  be  used  for  a single  point-design  estimate- 
using  manually  derived  inputs,  the  basic  costing  concept  requires  an  interface  with 
computerized  design  synthesis  programs  as  a primary  source  of  cost-related 
explanatory  variables. 

Section  2 describes  the  multistation  structural  synthesis  program  including  a 
description  of  input  requirements.  Section  3 describes  the  secondary  structure 
synthesis  program  together  with  inputs  required. 
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SECTION  II 

MULTISTATION  STRUCTURAL  SYNTHESIS  PROGRAM 

Structural  synthesis  is  a way  of  satisfying  the  design  problem  of  defining  a piece  of 
structure  that  fulfills  requirements  of  strength  geometry  and  other  criteria.  It 
combines  material  properties,  structural  analysis  techniques  and  loading  environ- 
ments to  produce  a consistent  design.  In  broad  terms  structural  synthesis  is  the 
replacement  of  traditional,  laborious  hand  calculations  with  an  automated  series  of 
logical  steps.  It  offers  real  advantages  in  solution,  speed,  and  accuracy.  Mathe- 
matical optimization  techniques  have  been  incorporated  that  are  untractable  by  hand 
calculations. 

A multi-station  approach  is  used.  Design  synthesis  proceeds  systematically  from 
root  to  tip,  in  discrete  steps,  usually  at  rib  locations,  in  a two  phase  system.  In  the 
first  phase  of  the  synthesis  process,  a set  of  initial  member  size  estimates  is  ana- 
lyzed. Margins  of  safety  are  computed.  Thickness  variables  of  all  elements  are  ad- 
justed by  iterative  steps  until  each  element  has  a zero  margin  of  safety  or  until  a 
minimum  gage  constraint  is  encountered.  The  second  phase  seeks  to  maximize 
margins  of  safety  by  refinement  of  element  geometry  while  holding  structural  weight 
constant.  When  this  has  been  accomplished  the  design  is  recycled  through  phase  one 
to  further  refine  structural  weight.  This  logic  is  repeated  until  satisfactory  conver- 
gence is  obtained. 


The  fact  that  margin  of  safety  maximization  rather  than  weight  minimization  is  used 
in  the  second  phase  permits  use  of  unconstrained  function  optimization  methods. 
Major  advantages  of  this  approach  are:  the  member  sizes  can  be  constrained  to  lay 
within  practical  limits  of  material  sizes  and  manufacturing  capability;  multiple  fail- 
ure modes  may  be  taken  into  consideration  for  each  structural  element;  positive 
margins  of  safety  are  always  maintained  so  that  a satisfactory  design  - from  the 
strength  point  of  view  if  not  the  weight  - is  available  at  the  completion  of  each 
iteration. 

An  accurate  representation  of  geometry  is  permitted  by  defining  discrete  nodes  on 
the  contour  of  the  surface.  The  calculation  of  internal  loads  distribution  is  im- 
proved over  previous  programs  by  incorporation  of  methodology  for  analysis  of  a 
multi-cell  box  beam.  Complex  ^ending,  shear  and  torsional  loads  may  be  applied. 
Axial  loads  and  shear  flows  are  computed  for  each  node  point  and  panel.  Beams  are 
limited  to  a maximum  of  four  cells.  The  discrete  nodes  used  in  defining  the  contours 
are  also  used  as  elemental  centers  of  mass.  A spinoff  of  this  modeling  scheme  is 
the  ability  to  represent  surfaces  using  the  dated  constructional  mode  of  concentrating 
the  bending  material  in  the  spar  caps. 
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Figure  3.  Flow  Diagram  for  Synthesis  Program  Inputs 
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Figure  4.  APAS  Function  Flow  Chart. 


structure  is  placed  inside  the  aerodynamic  envelope  at  each  control  station.  Member 
configuration  must  be  considered  and  selection  of  the  form  of  several  detail  com- 
ponents must  be  made.  Initial  estimates  of  member  sizes  are  satisfactory  because 
the  program  optimization  scheme  rapidly  converges.  A maximum  of  four  torsional 
cells  is  permitted  (Figure  6). 
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Figure  *5.  Box  Configurations . 

Section  contours  are  specified  by  locating  discrete  nodal  points  on  the  surface 
igure  1).  Panels  connect  these  nodes  and  form  the  box  covers  and  spar  shear 
v.ebs.  Each  panel  requires  a selection  of  panel  types,  initial  size  estimates  and  a 
‘-i/.e  range  that  bounds  the  tvpe  dependent  design  variables  The  size  range  may  be 
used  to  introduce  various  design  constraints.  The  lower  limit  or  thickness  may  be 
•••el  by  machinining  minimums  or  handling  requirements  while  the  upper  limit  may  be 
set  by  available  extrusion  sizes.  Another  practical  aspect  of  the  constraint  feature 
. the  ability  to  hold  certain  variables  constant.  A uniform  stringer  pitch  can  be 
maintained  or  stringers  kept  at  a constant  height. 

sample  panel  specification  is  shown  in  Figure  8 and  is  typical  of  the  type  of 
decisions  requi red  to  define  panel  components.  Spar  caps,  spar  webs,  rib  caps 
and  rib  webs  require  a similar  type  of  specification. 
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Figure  7.  Typical  Section  Definition 
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Rib  spacing  is  not  an  optimized  variable.  Practical  experience  has  been  that  rib 
locations  are  almost  universally  dictated  by  design  requi rements , usually  the  neces- 
sity of  providing  support  structure.  Pib  spacing  may  be  set  at  each  control  station 
and  is  held  constant  until  reset  at  another  control  station.  Spar  locations  must 
also  be  re-selected.  In  general,  exploration  of  the  effects  of  geometric  variations 
of  the  major  structural  items,  ribs  and  spars,  requires  multiple  runs. 


2.2  EXTERNAL  LOADS 


External  forces  for  up  to  six  loading  conditions  may  be  considered.  Loads  may  be 
lescribed  at  a maximum  of  20  stations  and  may  contain  three  shears  and  three 
moments.  Two  load  sets  may  be  applied  at  particular  stations  to  represent  discrete 
load  changes  that  might  arise  from  concentrated  loads,  external  stores,  etc. 


Ihe  program  processes  each  load  condition  at  each  load  stai’-  'n.  Straight  line  inter- 
polation is  used  to  estimate  load  values  for  ribs  between  stations  of  load  application. 

‘ „ iidltlcij*  ItvaUwl  In  •*  roi«pfi'U*ly  fcniiitiTOf  ftthfan.  I’rrU  of  il*.  feat 
sensitive  to  a particular  loading  condition  are  designed  to  that  condition  rather  than 
to  an  envelope. 


Special  loads  data  is  required  for  Right  safety  analysis.  Fail  safe  analysis  requires 
a scaled  down  load  factor.  Cyclic  loads  information  is  required  to  investigate  safe 
life  and  fatigue  effects. 


2.2  INTERNAL  LOADS  DISTRIBUTION 


It  is  possible  to  compute  a reasonable  estimate  of  the  internal  loads  distribution 
for  many  aircraft  components  by  simp  e beam  theory.  The  structure  under  consi- 
deration mist  be  capable  of  being  mrE  led  as  a beam.  It  should  have  continuity  and 
;t  wtwmnaWy  weU  Wlr#o  elastic  axis  Tbb  swrtw  landing 

md  torsional  moments  of  inertia  and  cross-section  areas  are  computed  in  the 
program.  External  loads  are  then  applied  at  each  station,  transierred  to  the  section 
elastic  axis  and  used  to  derive  a set  oi  internal  loads  for  each  external  loading 
cond  tion. 


Distribution  of  shears  around  the  section  is  found  by  algebraic  summation  of  the 
circular  auear  produced  by  the  applied  torsional  moment  (T/2A)  and  "he  shear  pro- 
duced by  direct  loads.  Direct  shear  is  distributed  according  to  YQ/I  relationships. 
The  effect  of  taper  in  relieving  direct  shear  stress  is  recognized  by  the  program. 
Shear  Row  varies  across  the  panel  elements;  the  shear  flow  magnitude  at  the  center 
of  the  panel  is  used  as  representative  for  analysis  purposes.  (Panel  elements  are 
the  member  that  connect  the  nodes  and  supply  the  structural  material.) 


y 


-« 


Stresses  in  the  longitudinal  direction  results  from  the  application  of  bending 
moments  and  axial  loads.  The  distribution  of  axial  loads  is  proportional  to  the 
cross-sectional  area.  Stresses  due  to  the  application  of  bending  moments  are 
calculated  under  the  assumption  that  plane  areas  remain  plane  after  bending 
(Mc/I).  Computations  are  made  at  each  node.  Stress  level  used  for  element 
analysis  is  the  average  value  of  stress  for  the  nodes  on  either  side  of  the  panel. 


A conservative  estimate  of  rib  load  is  made  by  assuming  that  the  change  in  shear 
between  stations  is  attributable  to  shear  load  applied  to  the  intervening  rib.  This 
shear  force  is  applied  to  the  rib  and  reacted  with  a VQ/l  shear  flow  distribution. 


The  distribution  of  material  leads  to  the  stiffness  characteristics  of  the  structure 
and  to  the  distribution  of  internal  loads.  As  the  design  is  refined  by  the  optimi- 
zation process,  changes  occur  in  the  amount  and  arrangement  of  the  structural 
material.  Consequently  internal  loads  analysis  is  a requi:ed  step  in  each  design 
iteration. 


2.4  MATERIALS 


A library  of  materials  is  available  to  the  program.  Mechanical  properties  are 
given  as  a function  of  temperature  for  the  aluminum  alloys  2219-T  87 , 2024-T6, 
2024-T851,  7075-T651;  titanium  alloys  Ti-8A1  -IMo-IV,  Ti-GA1-4V;  Rene  41  and 
an  arbitrary  boron  epoxy  composite  material  (40%  ± 45",  50*7  0C,  10%  90°).  An 
additional  array  is  left  vacant  for  temporary  storage  of  a material  not  provided  in 
the  library.  A linear  interpolation  scheme  is  used  to  estimate  mechanical  property 
values  at  the  structural  temperature  of  interest.  A similar,  although  smaller, 
library  of  mechanical  properties  of  honeycomb  core  is  also  provided.  Aluminum 
ceres  range  in  density  from  2.  6 Ib/cu.  ft.  to  8. 1 Ib/cu  ft.  ; steel  cores  from 

4.5  lb/cu.  ft.  to  10.7  !b/cu.  ft. 


Material  selection  requires  . ns i deration  of  many  design  parameters.  Examples 
of  some  parameters  are  ultimate  strength,  fatigue  resistance,  and  fracture 
characteristics.  Candidate  materials  may  be  compared  by  using  a multiple  run 
strategy. 


2.5  COMPONENT  ANALYSIS 


This  section  describes  the  primary  structural  box  components  that  are  integrated  to 
form  a load  carrying  structure.  Covers  basically  provide  material  to  resist  bending, 
spars  to  resist  shear  and  ribs  to  stabilize  the  section  and  provide  points  for  intro- 
duction of  concentrated  loads  into  the  structural  system.  Each  of  these  components 
may  be  fabricated  by  several  construction  methods.  Method  selection  may  be  a 
function  of  geometric  and  mission  parameters. 
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2.5.1  COVERS.  Ihree  general  structural  concepts  lor  cover  construction  are 
considered: 


Plate  Stringer  - Plate  stringer  construction  is  found  on  many  aerospace  structures. 
Skin  panels  are  supported  by  closely  spaced  stiffening  elements.  These  stiffeners 
may  be  a variety  of  shapes  and  may  be  Integral  with  the  skin  or  attached  by 
mechanical  means.  Popular  stringer  shapes  are  blade,  zee,  tee,  jay  and  hat. 
Ultimate  strength  of  the  composite  plate  stringer  elements  under  combined  com- 
pret  sive  and  shear  loads  is  the  usual  criteria  for  margin  of  safety  computation. 

Back  ed  skin  and  partially  failed  elements  are  not  precluded  by  the  analysis  methodo- 
logy. Stresses  are  also  compared  with  ultimate  material  properties. 


Multi  -web  - Multi -web  construction  is  often  used  on  nigh  load  factor  aircraft  such 
as  military  combat  types.  It  permits  concentration  of  the  bending  resistant  material 
as  close  as  passible  to  the  external  contour.  The  high  loads  felt  by  this  type  of 
structure  lead  to  relatively  thick  skins. 


Deflections  that  alter  the  section  contour  are  aerodynamkally  undesirable  and 
buckles  in  relatively  heavy  plate  tend  to  precipitate  cracks,  therefore  initial 
buckling  in  combined  shear  and  compression  s considered  to  be  a mode  of  failure. 
Other  failure  modes  are  net  stress  in  tension  and  compression.  A minimum 
thickness  cutoff  is  also  tested. 


FuU  Depth  Sandwich  - This  method  of  construction  is  used  on  thin  sections.  Core 
weight  is  usually  prohibitive  if  maximum  section  thickness  is  on  the  order  of  4-5 
inches.  The  skin  is  assumed  to  be  fully  stabilized  and  critical  failure  modes  are 
based  on  net  stresses  in  tension  and  compression  with  shear  interaction.  A mini- 
mum skin  thickness  may  be  specified. 


2 5.2  SPAR  CAPS.  The  majority  of  modem  designs  use  a load  carrying  skin. 
Hence  the  significance  of  the  spar  cap  in  resisting  bending  loads  is  greatly  reduced. 
Design  requi rements  evolve  to  supporting  the  skin  in  crippling  modes  and  providing 
shear  continuity.  The  limits  of  ultimate  mechanical  properties  must  also  be 
observed. 


An  exception  to  the  above  was  a recently  icountered  hot  structure  where  the 
bending  material  was  concentrated  in  the  opar  caps  while  covers  were  relatively 
unconstrained  to  permit  thermal  growth.  This  case  can  be  treated  by  the  program. 
Caps  compatible  with  built  up,  integral,  truss  and  corrugated  web  type  spar  are 
analyzed. 


The  basic  structural  shape  is  a T",  for  which  minimum  and  maximum  thickness 
and  leg  lengths  must  be  input.  (An  angle  can  be  represented  by  inputing  zeros  for 
thickness  and  length  of  an  outstanding  leg.  i Thickness  of  the  attaching  web  or  skin 
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is  also  required  (Figure  9).  The  size  range  is  usually  a function  of  manufacturing 
practicalities.  Available  mill  sizes,  minimum  machinable  gauge,  minimum  space 
for  fastener  attachment  and  handling  requirements  are  typical  examples  of  external 
constraints. 


FWD 


Figure  9.  Typical  Spar  Cap  Model  - Built  I'p  Web. 

Failure  modes  considered  are  gross  stresses  in  tension  or  compression , flange 
crippling  and  flange  net  shear.  Column  mode  of  failure  must  also  be  considered  for 
the  cap  when  truss  type  spars  arc  used,  it  is  conservatively  assumed  that  foi  thi.- 
case  column  length  is  equal  to  the  rib  pitch.  .Symmetrical  sections  must  be  user'.. 
Loads  are  transferred  to  the  cips  a;  discrete  points,  rather  than  continuously  t.-> 
with  other  designs;  therefore,  the  !ocal  loading  condition  is  a combination  axial  load 
and  shear.  Local  shear  flow  is  estimated  to  be  3 times  the  average  shear  How  for 
the  bay.  Either  built  up  or  integral  truss  type  spars  may  be  considered. 

2.  5.3  SPAR  WEBS.  Spar  webs  are  the  direct  shear  carrying  components  of 
structural  boxes.  They  also  contribute  to  the  support  of  torsionally  induced  shears. 
The  shear  carrying  capability  may  also  be  furnished  by  a diagonal  member  carying 
tension  or  compression.  The  method  of  analysis  arid  applicable  failure  modes  are 
greatly  influenced  by  the  type  of  construction.  Four  distinct  spar  web  types  have 
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been  identified  and  provided  for. 


Built  l p Web  Type  - The  spar  web  is  assumed  to  have  stiffeners  on  one  side  of  the 
web.  AH  stiffener  legs  are  of  equal  thickness.  Any  cross  section  of  four  or  less 
rectangular  elements  may  be  used  (Figure  10).  Web  and  stiffener  must  be  of  the 
same  material  Presently  provisions  for  2024-T3  and  7075-T6  aluminum  alloys 
are  included.  The  method  of  analysis  used  is  that  given  in  Reference  2.  Equations 
have  been  fitted  to  all  pertinent  curves.  Spar  web  allowable  stress  is  a function  of 
the  support  offered  by  the  surror  .ling  structure  and  the  material  properties. 
Stiffeners  are  checked  for  column  failure  and  forced  crippling. 


Integral  Web  - Ihe  integral  stiffener  analysis  assumes  a rectangular  stiffener 
cross  section  with  stiffener  thickness  equal  to  web  thickness  (I'igure  11).  An 
analysis  routine  was  established  based  on  the  methods  and  curves  of  Reference  3. 
Best  results  are  obtained  with  materials  for  which  the  ratio  of  F /F  is  approxi- 
mately 1.  5.  The  allowable  ultimate  web  stress  may  be  either  bucideiUor  unbuckled 
depending  upon  the  geometric  and  material  property  parameters. 


Corrugated  Web  - The  corrugated  web  consists  of  a series  of  circular  segments 
joined  in  a continuous  fashion  to  rigid  caps  (Figure  12).  There  are  no  theoretical 
material  restrictions,  however,  it  should  be  noted  that  the  most  practical  methods 
of  construction  require  some  type  of  fusion  process.  Usually  a weldable  material 
is  considered.  Analysis  methodology  is  given  in  Reference  4.  Corrugated  con- 
<A,'  iction  is  considered  to  have  Railed  when  buckling  either  local  or  general  occur. 
In  addition,  net  shear  stress  is.  checked. 


Truss  - A truss  arrangement  cun  serve  the  s:une  functional  shear  carrying  role  as 
a shear  web.  Application  takes  the  form  of  diagonal  axial  members  carrying  either 
tension  or  compression  (Figure  13).  It  is  often  difficu’t  to  insure  that  only  tension 
will  be  felt  bv  a particular  member  since  at  least  pai  dal  load  reversability  is  very 
common.  Each  member  is  analyzed  for  the  maximum  load  in  compression. 

Diagonal  cross  section  may  be  cruciform,  square  solid,  round  solid,  or  tube.  A 
representative  truss  angle  may  be  input  or  a default  value,  one  radian,  used. 
Failure  modes  are  column  buckling  for  all  shapes  and  also  local  crippling  for  the 
non-compact  sections.  'Ihis  analysis  is  applicable  to  both  integral  and  built  up 
truss  webs. 


2 5 4 RIBS'.  The  functional  role  of  rib  element:,  is  dependent  upon  the  structural 
concept  being  considered.  In  general  ribs  serve  to  help  maintain  section  contour 
support  air  loads  and  redistribute  shears  and  concentrated  loads.  The  most  signi- 
ficant factors  on  rib  design  are  usually  the  redistribution  of  shears  and  the  impact 
of  concentrated  loads. 


Formed  sheet  metal  ribs,  corrugated  ribs,  ribs  with  built  up  or  integral  webs  and 
built  up  or  integral  truss  ribs  are  the  structural  types  available  in  the  rib  library 
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Figure  12.  Corrugated  Web 
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Figure  14).  Rib  web  shears  are  determined  by  a VC-i/l  analysis  where  V is  the 
span  wise  change  in  shear  for  one  rib  bay.  The  rib  web  is  checked  for  the  maxi- 
mum shear  flow  at  the  station.  When  full  depth  sandwich  construction  is  used  the 
core  material  serves  as  a form  ol  rib. 

Formed  Sheet  Metal  Rib  - The  formed  sheet  metal  rib  is  assumed  to  be  constant 
thickness  including  the  rib  flange  or  "cap”.  Failure  modes  for  this  rib  are  shear 
buckling,  net  shear  stress  exceeding  material  properties,  crushing  due  to  bending 
and  shear  in  the  flange. 

Corrugated  Rib  - Corrugated  ribs  are  formed  by  a continuous  wave  of  circular 
segments.  Web  thickness  is  uniform.  Rib  caps  offer  continuous  support.  Net 
shear,  shear-compression  interactive  buckling,  minimum  gage,  and  net  shear 
through  the  flange  are  the  modes  of  l';<ilure  considered. 

Built  Up  Web  - The  built  up  rib  is  very  similar  in  construction  to  the  built  up  spar 
Analyses  methods  are  also  similar.  Allowable  stress  in  the  web  is  a function  of  the 
support  provided  by  the  stiffeners  and  caps.  A typical  stiffener  is  checked  for 
column  failure  and  forced  crippling. 

Integral  Web  - The  integral  web  rib  is  virtually  identical  to  the  integral  web  spar. 
Failure  mode  may  be  either  buckled  or  unbuckled  depending  upon  the  fixity  provided 
by  the  surrounding  structure.  The  same  analysis  routines  are  used  for  both  spar 
and  rib  webs. 

Truss  - Truss  rib  elements  may  be  selected  from  a library  of  four  simple  shapes 
Buckling  and  crippling  failure  modes  are  investigated  for  a typical  element.  A 
representative  truss  angle  may  be  input  or  a default  value  ol  one  radian  will  be  used. 
Integral  and  built  up  truss  ribs  are  treated  identically  by  the  program. 

Core  - Core  material  is  available  in  discrete  increments  rather  than  the  continuum 
of  gages  available  to  the  other  modes  of  construction.  Analysis  of  the  complex 
failure  modes  of  sandwich  core  are  loregone  and  load  magnitude  is  used  to  select  im- 
proper core  density. 

2.6  FLIGIH  SAFETY  ANALYSIS 

Fail  safe,  sate  life  and  fatigue  analysis  are  different  methods  of  approaching  the 
problem  of  structural  reliability  for  the  expected  life  and  usage  of  the  vehicle.  Th 
fail  safe  concept  relies  on  redundancy  of  structure  to  support  the  required  load 
after  a critical  member  has  failed.  Presumably  the  vehicle  will  be  able  to  complete 
its  flight  and  the  damage  will  be  discovered  and  repaired.  Iitmage  could  be  the 
result  of  fatigue,  or  combat.  Fail  safe  philosophy  cannot  be  applied  to  all  de' !gns 
because  redundant  construction  is  not  always  feasible.  Fatigue  analysis  estimates 
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the  life  of  the  structure  under  environmental  conditions  of  temperature,  local 
geometry  (stress  concentrations)  and  cyclic  loading.  The  usual  criteria  examines 
the  number  of  cycles  or  structural  life  until  a crack  appears.  A fair  amount  of  data 
in  the  form  of  S-N  curves  has  been  established  on  fatigue  properties  of  various 
metallics.  The  safe  life  concept  recognizes  that  flaws  are  inherent  in  virtually  all 
metallic  structures  as  a result  of  the  manufacturing  process.  Flaws  down  to  a 
certain  size  can  be  eliminated  by  inspection  procedures  or  proof  testing.  The 
question  is,  what  is  the  life  of  the  structure,  given  an  initial  flaw  size  and  the 
design  spectrum  of  loads.  Critical  stress  intensity  factor  (K  ) and  crack  growth 
rate  (da/dn)  have  been  found  to  be  significant  parameters  in  cfetermining  safe  life. 


A scatter  factor  of  from  2 to  4 is  usually  required  on  safe  life  and  fatigue  analysis. 
Safe  life  and  fatigue  analyses  while  serving  the  same  general  purpose  are  not 
mutually  exclusive  and  are  sometimes  both  requested  by  service  agencies.  For 
operation  of  this  program  it  is  anticipated  that  the  user  will  specify  the  applicable 
type  ot  analysis  and  provide  the  required  inputs.  Program  output  will  be  in  the 
form  of  fail  safe  margins  of  safety  for  fail  safe  checks;  and  service  or  design  life 
for  fatigue  and  safe  life  analysis. 


A reasonable  indication  of  life  and  damage  resistance  can  be  obtained  if  the  selected 
analysis  is  applied  at  two  or  three  critical  location s on  the  structure.  For  a wing  or 
other  lifting  surface  the  tension  skin  at  the  root  and  one  or  two  intermediate  stations 
toward  the  tip  are  likely  points  of  investigation.  A pressurized  structure  should  be 
checked  at  a location  where  bending  and  pressurization  stress  add. 


2.6.1  SAFE  LIFE.  Safe  life  design  requires  sufficiently  low  design  stresses  that 


catastrophic  failures  of  critical  structural  components  will  not  occur  during  a 
specified  service  life  due  to  the  growth  of  flaws  or  defects  pre-existent  in  the 
structure.  The  concept  requires  that  the  structure  still  support  an  appreciable 
even  if  fairly  long  cracks  are  pie.sent  and  that  cracks  will  be  detected  with  100 
percent  certainty  before  they  can  grow  to  a dangerous  length. 


oad 


Plain  strain  fracture  toughness,  K is  a key  parameter  in  evaluating  resistance  to 
catastrophic  facture.  The  stress  intensity  factor,  K is  compared  to  K after  com- 
pletion of  each  block  of  cycles  (representing  on  flight).  Crack  growth  can  be  com- 


pared with  inspection  criteria. 


Safe  life  is  the  life  for  initial  defects  to  grow  to  a critical  size  for  catastrophic 
failure.  Typically  the  design  requirements  will  be  input  as  conditions  - usually 
points  of  the  flight  envelope  representing  the  corners  of  the  V-N  diagram  plus  a 
ground  condition  and  possibly  a few  maneuver  conditions. 


As  an  illustration  assume  the  following  design  conditions  are  applicable: 
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Condition  1 

Maximum  pull-up  (Nz  max.) 

2 

Buffet 

3 

Pushover 

4 

Rolling  pullout  (assymetric) 

5 

Taxi 

Assume  each  flight  may  be  characterized  as: 

20  cycles 

100%  of  Cond.  1 

40  cycles 

l/Nz  max.  Cond.  1 

15  cycles 

100%  of  Cond.  2 

10  cycles 

1.5/Nz  max.  Cond.  1 

10  cycles 

50%  Cond.  4 

2 cycles 

50%  Cond.  3 

20  cycles 

100%  Cond.  5 

(This  information  will  provide  stress  levels.) 

A desired  number  of  flights  is  also  required  and  usually  one  of  the  specifications 
is  1000  flights.  Initial  flaw  size  will  be  an  input  and  a function  of  the  inspection  or 
testing  program.  The  flaw  growth  and  stress  intensity  factor  will  be  computed  at  the 
end  of  each  flight.  The  allowable  number  of  flights  will  be  determined  when  the  K 

Q 

is  equal  to  K critical  - or  alternatively  if  the  ,<■  ified  number  of  flights  is  reached 
before  exceeds  critical  the  design  can  be  c msidered  satisfactory. 

One  station  at  a time  will  be  analyzed  and  the  results  will  be  informational 

desig  is  , .cceptable;  or  N flights  can  be  completed  before  the  stress  intensity  factor 
exceeds  a critical  value. 

2.6.2  FAIL  SA  TE . Fail  safe  design  requires  that  the  failure  of  any  single  structur- 
al component  will  not  degrade  the  strength  or  stiffness  of  the  remainder  of  the 
structure  to  the  extent  that  the  vehicle  cannot  complete  the  mission  at  a specified 
percentage  of  limit  loads. 

Application  of  this  philosophy  leads  to  redundant  designs  offering  many  possible 
load  paths.  At  preselected  critical  locations  members  are  assumed  to  have  failed 
and  internal  loads  in  the  local  structure  are  computed  for  fail  safe  applied  load. 

Fail  safe  margins  of  safety  may  also  be  computed. 

For  commercial  transport  aircraft  where  safety  is  of  utmost  concern,  fail  safe 
capability  is  provided  to  the  greatest  possible  extent.  For  military  aircraft  where 
performance  is  of  primary  concern,  fail  safe  capability  is  not  provided  where 
weight  penalties  would  result.  Other  methods  described  in  the  section,  safe  life 
and  fatigue  analysis,  are  relied  upon  to  estimate  flight  safety  in  the  case  of  military 
aircraft. 
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2.6.3  FATIGUE.  Aircraft  structural  components  loaded  in  tension  are  often 
designed  to  preclude  fatigue  failure.  Components  whose  design  allowables  and 
weights  are  influenced  by  fatigue  are  lower  wing  skins  and  pressurized  fuselages. 

Fatigue  analysis  is  a complex  procedure  that  requires  substantial  amoi  ts  of  data 
and  criteria.  It  is  generally  performed  at  critical  locations  on  the  components  of 
interest.  To  facilitate  synthesis  sr  ,t  .al  simplifying  assumptions  based  on  design 
experience  can  be  made  to  reduce  he  data  required  and  the  work  involved  to  provide 
fatigue  criteria.  The  fatigue  requirements  can  be  provided  in  summary  form  offering 
the  comparison  of  design  service  life  and  stress  level. 

The  methodology  for  development  of  curves  portraying  Ft/Ftu  vs  design  sendee  life 
is  given  in  Reference  5.  Results  for  several  cases  of  interest  are  also  presented. 
These  cases  are  curved  fitted  and  incorporated  in  a fatigue  subroutine. 

Fatigue  life  is  the  life  an  unflawed  structural  component  to  the  initiation  of  visible 
fatigue  cracks.  It  is  synonymous  with  service  life. 

2.7  AEROELASTIC  EFFECTS 

Aeroelastic  effects  on  the  loads  distribution  have  been  included  in  previously 
prepared  Convair  programs  (References  6 and  7)  Use  of  these  programs  is 
recommended  to  generate  the  external  load  data.  Other  aeroelastic  effects  of 
interest  are  straight  wing  divergence  and  flutter  and  aileron  reversal.  The  steps  of 
the  wing  sizing  procedure  generate  mass  and  stiffness  data.  This  data  is  required 
to  prepare  the  math  model  that  is  required  for  the  aeroelasuc  analysis. 

Results  of  this  analysis  will  indicate  divergence  and  aileron  reversal  speeds  and 
frequencies  of  the  first  bending  and  torsional  modes.  These  results  may  be  com- 
pared w.-iF  design  criteria  and  previous  experience  with  similar  designs  and  a 
judgement  as  to  acceptability  may  be  made. 

If  the  design  is  unacceptable  corrective  action  is  required.  The  appropriate  action 
is  seldom  obvious.  Many  possible  measures  may  offer  solutions.  Some  possibilities 
are  criteria  revision,  addition  of  damping'  systems,  revision  of  structural  arrange- 
ment, material  change  and  redistribution  of  material.  The  weight  and  cost  impact 
of  these  various  approaches  and  the  methodo’ogy  to  achieve  them  suggest  complex 
questions  that  are  best  considered  outside  the  synthesis  loop. 


Presently  aeroelastic  effects  are  studied  in  an  external  program.  Driving  inputs  are 
supplied  by  synthesis  and  results  are  evaluated  by  disciplinary  specialists. 


2.8  WEIGHTS 

Since  all  of  the  geometric  information  is  availab':„  to  the  program,  computation  of 
volume  and  weight  of  the  structural  material,  is  a relatively  simple  step.  This 
computed  weight  will  account  for  only  a fraction  (historically  approximately  l/2) 
of  the  actual  finished  weight.  Items  will  exist  in  every  structure  that  ire  insensitive 
to  load  or  are  impossible  to  identify  for  synthesis.  Weight  correlation  factors  are 
required  to  adjust  the  computed  weight  to  a more  represented  value. 

2.9  A PAS  PROGRAM  INPUTS 

Hie  multi -station  structural  synthesis  program  performs  analysis  and  redesign 
operations  on  one  station  at  a time  in  a systematic  fashion  from  root  to  tip.  Each 
loading  condition  is  processed  and  the  full  complement  of  structural  elements  at 
that  station  are  satisfactorily  optimized  before  subsequent  stations  are  considered. 
Input  requirements  for  program  execution  are  identified  by  the  set  of  input  caids 
described  below. 


CARD  No.  1 
Format  (413) 


OUTPUT  CONTROL  CARD 


KOU  T1 , KOU  T2 , KOU  T3 , KOUT4 

Allows  user  to  output  various  parts  of  the  output  to  different  output  devices.  A 
blank  card  defaults  all  output  to  device  6 (printer). 


CARD  No.  2 


ITERATION  AND  TOLERANCE  CONTROL  CARD 


Format  (4  (15,  5X) , 4F10. ) 

ITl.I  F2 , IT3 , IT4,  EPS1,  EPS2 , EPS3 , EPS4 

Allows  user  to  set  '.units  on  accuracy  of  solution  and  maximum  passes  through 

various  iteration  loops.  A blank  card  defaults  all  terms  to  suggested  values. 

IT1  Iteration  count  limit  on  overall  redesign/optimization  procedure 

default  value  is  (5) . 

IT2  Iteration  count  limit  on  fully  stressed  redesign  process. 

IT3  Iteration  count  limit  on  Fletcher  Powell  optimization  procedure. 

IT4  Not  used. 

EPS1  Tolerance  on  redesign  margins  of  safety  for  each  redesign  cycle 
(Default  value  is  .001) 


EPS2  Tolerance  on  final  margins  on  safety.  At  least  one  non-minimum  gage 
element  of  each  symmetry  group  for  at  least  one  load  condition  will  have 
a margin  of  safety  which  satisfies 

MS  £ EPS2 

(Default  value  is  .01). 

EPS3  Tolerance  on  optimization  function  decrease  in  Fletcher  Powell 
minimization  technique. 

(Default  value  is  .001). 

EPS4  Toierance  on  design  variable  variation  in  Fletcher  Powell  minimization 
technique 

(Default  value  is  . 001). 


CARD  No.  3 


TITLE  CARD  (2  ca  ds) 


Format  (8A10) 
TITLE  (I)  1 = 1,  16 


CARD  No.  4 


CASE  CONTROL  CARD 


Format  (1013) 

KEY1 , KEY2,  KEY3,  KE  Y4,  KEY5,  KEY6,  KEY7,  KEY8.  KFY9.  KEY10 

KEY1  Allows  user  to  select  one  of  four  geometry  input  subroutines 

KEY1  = 1 specifies  a general  input  subroutine  "GINPT]" 

which  is  described  below. 

= 2 not  currently  available 

= 3 not  currently  available 

= 4 noi  currently  available 


KEY2  Allows  user  to  select  one  of  four  external  loads  definition  subrouL 

KEY2  = 1 specifies  a load  input  subroutine  "LOADIN'-  describe 
below  by  which  the  user  inputs  all  loads. 

KEY2  = 2 not  currently  available 

= 3 not  currently  available 

= 4 not  currently  available 


KEY3  =0  for  fuselage  type  structures 
= 1 for  wing  like  structure 
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KEY4  Allows  user  to  select  how  many  locations  along  the  structure  he  wishes 
to  synthesize. 

KEY4  = 1 synthesize  at  every  rib/frame 

= 2 synthesize  at  every  other  rib/frame 
= 3 synthesize  at  every  third  rib/frame 


KEY5  Not  used 

KEY6  Allows  user  to  check  input 


KEY7 


KEY6  0 normal  mode 

1 check  input  and  quit. 

Rib/Frame  Type  Specification 
Rib  (KEY3  = 1) 


KEY7  = 1 corrugated  web 

= 2 integral  web 

- 3 built  up  web 

= 4 integral  truss 

= 5 built  up  truss 

Frame  (KEY3  = 0)  not  currently  available 

KEY8  Unsupported  panel  width  in  terms  of  percent  of  panel  element  width. 
Used  for  panel  types  10,  11  and  12. 

KEY9  Not  used 

KEY10  Not  used 


CARD  No.  5 


MATERIAL  SPECIFICAT  ION  CARD 


Format  (15,  5x,  F10.0) 


IMAT,  FTEN 

ima  r 


1 user  supplied  metallic  material 

2 AL-2219-T87 

3 TI-8AL-1M0-1V 

4 AL-2024-T6 

5 AL-7075-T6 

6 Inconel  718 

7 Inconel  625 

8 TI-6AL-4V 

9 AL-2024-T851 
10  Rene'  41 


FTEN 


CARD  No.  5A 


11  User  supplied  composite  material 

12  NARMCO  5505  Boron-Epoxy 

13  N \RMCO  5206  Graphite -Epoxy 

Factor  applied  to  ultimate  tension  allowable  to 
account  for  fatigue  fracture  mechanics,  etc. 

FT  = FTEN  * FTU 

USER  SUPPLIED  METALLIC  MATERIAL 
PROPERTIES.  (USED  WHEN  IMAT  = 1) 


This  option  allows  the  user  to  input  a material  not  found  in  the 
material  table  above.  The  properties  may  be  input  as  a function 
of  temperature  at  nine  specified  temperatures. 

5A1  Materia]  title  card 
Format  fJA10) 

5A2  NTEMP 

Format  (110) 

NTEMP  is  the  number  temperature  points  specified. 

5A3  Room  temperature  properties  cards  (2  cards) 

Format  (3F10.0,  2F20. 0/F20.  0, 3F10. 0) 

FTU  FCY  FSU  EC  E G RHO  F07  EN 

5A4  Temperature  factors  cards.  (9  cards  maximum) 

Format  10 F8. 0) 

TEMP,  FFTU , FFCN,  FFSU , FEC,  FE,  FG,  FRHO, 
FF07 , FEN 

CARD  No.  5B  USER  SUPPLIED  COMPOSITE  MATERIAL  PROPERTIES 
(Used  when  IMAT=  11) 

This  option  allows  the  user  to  input  a composite  material 
not  found  in  the  material  table. 

5B1  Material  title  card 
Format  (8A10) 

5B2  Room  temperature  properties  cards  (2  cards) 

Format  (3F20.0,  F10.0/6F10.  0) 

Ell,  E22,  G12,  U12,  DEN,  EPSAL1, 2,  3 , 4,  5, 

CARD  No.  6 GEOMETRY  CONTROL  CARD 
Format  (413) 

NODES,  NW EB , NLONG , NSTAG 


1 


sp  f 


Vin 


NODES 


Number  of  node  points  in  a cross-section  (must  be  at  least 
2 and  may  be  maximum  of  20) 


NWEB 


Number  of  interior  webs  (Max.  of  3) 


NLONG  Number  of  spar  caps/'ongerons  (Max.  of  10)  (must  be 
located  at  nodal  points) 


NSTAG  Number  of  geometry  control  stations  (Max.  of  20) 


CARD  No.  7 


GEOMETRY  CONTROL  STATION  HEADER  CARD 


Format  (6F10. 0) 


STAG,  FRSP.XLDRF,  ZLDRF,  GTRX.GTRZ 


STAG  Station  number  of  control  station 


FRSP  Rib/ Frame  Spacing 


XLDRF  Coordinates  of  input  loads 


7LDRF  Reference  Point 


GTRX  Taper  ratio  of  structure 


GTRZ  In  the  X and  Z directions 


CARD  No.  8 GEOMETRY  NODAL  COORDINATE  DECK 


Format  (15,  fX,2F10.0)  (1  card  for  each  node) 


NODE , X,  Z 


NODE  Node  number  start  with  1 and  number  consecutively 
around  the  structure  clockwise. 


Coordinates  of  node  point 


Repeat  Cards  7 and  8 for  each  geometry  control  station. 


CARD  No.  9 INTERNAL  WEB  LOCATION  CARD 


Format  (613)  (Omit  if  NWEB  is  0) 


IW7  , LW1 , IW2,  LW2 , IW3 , LW3 . 


Node  number  of  first  node  to  which  an  interior  web  is  attached 
(Web  No.  1) 


LW1  Node  number  at  other  end  of  web  No.  1 


Same  as  above  for  next  interior  web. 
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Same  as  above  for  last  interior  web. 


CARD  No.  10  SPARC AP/LONGE RON  NODE  IDENTIFICATION 

DECK 


Format  (4  (13, 2X,  F10.  0,5X) 


LI , DL1 , L2,  DL2,  etc.,  four  sets  per  card  - one  set  for  each  longeron 


Node  number  of  spar  cap/longeron  i 
Orientation  angle  of  spar  cap/longeron 


CARD  No.  11  SYMMETRY  GROUP  CONTROL  CARD 


Format  (313) 


NSGP,  NSGW,  NSGL 


NSGP  Number  of  panel  symmetry  groups 


NSGW  Number  of  interior  web  symmetry  groups 


NSGL  Number  of  spa  reaps  (Longeron)  symmetry  groups. 


CARD  No  12  PANEL  SYMMETRY  GROUP  IDENTIFICATION 
DECK  (1  card  per  symmetry  group) 


Format  (2413) 


NMAX.NSd)  I-l . NMAX 


NMAX  Number  of  path,  s in  this  symmetry  group 


NS(1) 


Panel  number  of  first.  p;tnel  in  the  symmetry  group  going 
CW  around  the  cross  section. 


NS  (2) 


Panel  number  of  next  pane! 


CARD  No.  13  INTERIOR  WEB  SYMMETRY  GROUP  IDENTIFICATION 
DECK  (1  card  per  symmetry  group) 


Follow  the  same  procedure  as  for  panel  symmetry  group  deck. 


CARD  No.  14  LONGERON  iSPARCAP)  SYMMETRY 
GROUP  IDENTIFICATION  DECK. 


Follow  same  procedure  as  for  panel  symmetry  group  deck. 
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CARD  No.  15 


STRUCTURAL  ELEMENT  CONFIGURATION  SPECIFICATION 
DECK 


This  deck  defines  the  structural  configuration  desired  for  each  panel,  web, 
or  longeron  (spar  cap).  A deck  is  required  for  each  symmetry  group  defined, 
an  additional  deck  is  required  for  each  panel,  web  or  longeron  (spar  cap) 
which  does  not  belong  to  a symmetry  group.  The  order  of  the  decks  is  as 
follows: 


PANELS 


Symmetry  group  1 thru  NMAX 

Panels  which  are  not  members  of  symmetry  groups, 
starting  with  the  lowest  such  panel  number  and  going 
CW  around  the  structure. 


WEBS 


Follow  same  procedure  as  for  panels. 


LONGERONS  (Spar-Caps) 


Follow  same  procedure  as  for  panels. 


CARD  No.  15A  TYPE  IDENTIFICATION  CARD 


Format  (213) 


ITYP,  IDSET 


ITYP  Structural  configuration  type  number 


IDSET 


(Opu  >nal)  Each  deck  may  be  given  an  ID  number  which  can  be 
referred  on  other  type  identification  cards  which  have  identical 
structural  element  configuration  specification  decks  When 
referring  to  an  ID  number  specified  on  a pre  vious  deck,  cards 
15b  through  15f  are  omitted.  (Note:  Web  Id's  may  not  refer 
to  panel  ID's , etc . ) 


CARD  No.  15b  T SPECIFICATION  CARD 


Format  (4F10.  0) 


T1.T2.T3.T4 


Initial  values  for  T variables. 


CARD  15c  TMIN  SPECIFICATION  CARD 


Format  (4F10.  0) 


TMIN1 , TM1N2 , TMIN  3 , T MIN 4 


§ 


r 


m 


B 


I 


wi 


m 


I 


Minimum  values  for  T variables  on  card  15b.  This  card  allows  the 
user  to  specify  minimum  gages  for  the  T variables.  He  may  also  fix 
a T variable  at  its  initial  value  by  setting  its  TMiN  to  0.  (Note:  At 
least  one  T variable  must  have  a non-zero  TMIN. ) 


CARD  15d 


B VARIABLE  SPECIFICATION  CARD 


Format  (4F10.  0) 

B1,B2,B3,B4 

Initial  values  for  B variables 


CARD  15e 


BMIN  SPECIFICATION  CARD 


Format  (4F10,  0) 

BMIN1,  BMIN2,  BMEN3,  BMIN4 


Minimum  values  for  B variables.  The  user  may  fix  any  or  all 
B values  to  their  initial  values  by  setting  their  respective  BMIN 
values  to  zero. 


CARD  15f 


BMAX  SPECIFICATION  CARD 


Format  (4F10.0) 

BMAX1,  BMAX2,  BMAX3,  BMAX4 


Maximum  limits  for  B variables.  These  values  are  ignored  if 
cor  re  soon  ding  B MIN  values  are  zero. 


CARD  No.  16  INPUT  LOADS  CONTROL  CARD 


Format  (15,  5X,  F10.  0) 

NCOND,  FULT 

NCOND  Number  of  loading  conditions  (Maximum  of  6). 
FULT  Ultimate  factor  of  safety. 

CARD  No.  17  LOAD  CONDITION  HEADER  CARD 


f ormat  (4A1 0 ,15,  5X,  F10.  0) 

Condition  title,  NSTAY,  press  condition  title  40  characters 
NSTAY  Number  of  stations  at  which  the  load  case  is  defined 
PRESS  Internal  pressure  used  for  fuselage  structures. 
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CARD  17a 


LOAD  COMPONENT  FACTOR  CARD 


Format  (8 F10.  0/F10.  0)  (2  cards) 


These  factors  applied  to  cards  17b) 


FLIN 


Factor  applied  to  input  stations  to  produce  station  numbers 
(e.  g.  , input  stations  could  be  entered  in  percent  fuselage 
length  (wing  span)  and  then  FLIN  would  be  actual  fuselage 
length  (wing  span) . 


Factor  applied  to  all  input  load  components 


FXS 
FZS 
FT  OR 


Factor  applied  to  AXIAL 
Factor  applied  to  XSHEAR 
Factor  applied  to  ZSHEAR 
Factor  applied  to  TORSION 
Factor  applied  to  XMOMENTS 
Factor  applied  to  Z MOM  ENTS 


loads 

loads 

loads 

loads 


FTEMP  Factor  applied  to  TEMPERATURE 


Note: 


Factored  load  components  are  used  by  the  program  as 
limit  design  loads . 


If  cards  17a  are  blank,  all  factors  are  defaulted  to  1. 


CARD  17b 


LOAD  POINT  CARD 


Format  (8F19.  0) 


STA  , AX,  SX,  ZS,  TOR,  XMOM,  ZMOM , TEMP 


xs 

ZS 

TOR 

XMOM 

ZMOM 

TEMP 


Station  Number 
Station  axial  load 
Station  X Shear  force 
Z SHEAR 
Torque 

Bending  moment  about  X axis 
Bending  moment  about  Z axis 
Structural  temperature  degrees  I’arenheit 
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SECTION  III 


SECONDARY  STRUCTURE  SYNTHESIS 


This  section  presents  a technical  discussion  associated  with  the  initial  development 
stage  of  a computer  program  for  estimating  the  geometry,  weights,  part  definition 
and  primary  cost  drivers  for  aerodynamic  surface  leading  edge,  trailing  edge  and 
tip  components.  The  aerodynamic  surface  structural  box  is  not  included  as  part  of 
the  program  development.  However,  the  output  is  designed  to  complement  structural 
box  analysis  to  provide  complete  coverage  of  the  aerodynamic  surfaces. 


The  discussion  is  presented  in  two  sections.  The  first  deals  with  the  geometry  and 
analysis  of  the  leading  edge,  trailing  edge,  and  tip  components  and  the  second  deals 
with  the  component  part  definition.  The  operation  and  output  of  the  geometry, 
analysis  and  part  definition  program  subroutines  are  designed  for  use  in  the  con- 
ceptual and  preliminary  design  phases  of  aircraft  development. 


3.1  TIP,  LEADING  AND  TRAILING  EDGE  ANALYSIS 


The  leading  edge,  trailing  edge,  and  tip  synthesis  modules  provide  the  capability  to 
analyze  the  aerodynamic  surface  structural  components  that  are  not  considered  as 
part  of  the  structural  box.  The  leading  edge  is  defined  as  being  forward  of  the  front 
spar  and  includes  the  fixed  portion  of  the  leading  edge  and  the  leading  edge  high  lift 
devices  (slats).  The  trailing  edge  is  defined  as  being  aft  of  the  rear  spar  and 
includes  the  fixed  trailing  edge,  foreflaps,  flaps,  ailerons,  rudder,  elevator,  and 
spoilers.  The  tip  is  defined  as  that  structure  outboard  of  the  structural  box  tip 
closing  rib. 


The  synthesis  includes  a definition  of  part  geometry  and  a detailed  stress  analysis 
that  determines  gages,  accounts  for  material  types,  and  sets  minimum  gage  con- 
straints. The  geometry  routines  provide  dimensional  input  to  the  stress  analysis 
routines.  The  geometry  and  stress  routines  output  includes  part  size  and  weight, 
as  well  as  parameters  for  the  part  definition.  A generalized  flow  of  the  leading 
edge,  trailing  edge,  and  tip  subprogram  is  shown  in  Figure  15. 


Ihe  analysis  utilizes  eight  geometry  routines,  three  stress  analysis  routines,  six 
supporting  routines,  and  two  calling  routines.  The  geometry  routines  are  for  flaps 
aileron,  rudder,  elevator,  slat  location,  slats,  fixed  leading  edge,  and  spoilers. 


The  stress  analysis  routines  include  foreflap,  spoiler,  and  one  which  analyzes  the 
flaps,  ailerons,  slats,  rudder,  and  elevator.  The  supporting  routines  derive  dimen- 
sions, material  properties,  and  general  analysis.  A discussion  of  these  routines  is 
included  in  the  following  paragraphs. 


* 
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SUPPORTING  SUBROUTINES: 
THICK  DENS 


Figure  15.  Leading  Edge  and  Trailing  Edge  Synthesis  Routines. 
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The  flap  geometry  routine  provides  flap  planform  dimensions  and  locations  from 
input  data.  The  flap  types  considered  are  simple  flaps,  and  single-slotted  and 
double-slotted  flaps.  In  the  case  of  single  or  double  slotted  flaps  the  ioreflap 
dimensions  are  computed  in  addition  to  the  main  flap  dimensions.  The  driving  para- 
meters in  determining  flap  dimensions  are  the  flap  area  to  wing  area  ratio,  flap 
chord  to  wing  chord  ratio,  and  flap  inboard  chord.  If  the  area  ratio  is  input  the 
flap  length  will  be  set  to  give  required  flap  area  The  flap  length  will  be  truncated 
at  the  wing  tip  or  the  inboard  edge  of  the  aileron.  The  flap  chord  is  set  by  the  ratio 
of  flap  chord  to  wing  chord.  If  the  ratio  is  zero  the  chord  is  assumed  to  be  85%  of 
the  distance  aft  of  the  rear  spar.  If  the  flap  chord  is  input,  the  value  of  flap  chord 
to  wing  chord  ratio  will  be  computed  for  use  in  determining  flap  dimensions.  The 
inboard  edge  of  the  flap  is  located  at  the  side  of  the  fuselage.  Flap  geometry  output 
consists  of  inboard  and  outboard  chords,  span  stations  of  the  flap  inboard  and  out- 
board edges,  and  the  flap  length. 


The  aileron  geometry  routine  provides  ai'eron  planform  dimensions  and  locations 
from  input  data.  The  outboard  edge  of  the  aileron  is  assumed  to  be  at  the  wing  tip 
and  the  inboard  edge  is  truncated  at  the  side  of  the  body  if  the  inboard  location  is  not 
specified  The  aileron  chord  is  computed  as  10%  greater  than  the  trailing  edge 
length.  If  the  inboard  edge  location  of  the  aileron  is  input  the  length  will  be  set  to 
provide  the  required  aileron  area.  Aileron  geometry  output  consists  of  inboard  and 
outboard  chords,  span  stations  of  the  aileron  inboard  and  outboarc.  edges,  and  the 
aileron  length. 


The  rudder  geometry  routine  provides  rudder  planform  dimensions  and  locations 
from  input  data.  The  rudder  extends  from  the  body  to  the  vertical  stabilizer  tip. 

The  rudder  chord  is  set  as  90%  of  the  distance  aft  of  the  vertical  stabilizer  rear 
spar.  Rudder  geometry  output  consists  of  inboard  ;md  outboard  chords,  span  stations 
of  the  rudder  inboard  and  outboard  edges,  and  the  rudder  length. 


The  elevator  geometry  routine  provides  elevator  planform  dimensions  and  locations 
from  input  data.  The  elevator  extends  from  the  body  to  the  horizontal  stabilizer  tip. 
The  elevator  chord  is  set  as  90%  of  the  distance  aft  of  the  horizontal  stabilier  rear 
spar.  Elevator  geometry  output  consists  of  inboard  and  outboard  chords,  span 
stations  of  the  elevator  inboard  and  outboard  edges,  and  the  elevator  length. 


The  slat  geometry  routine  comprises  two  separate  operations.  The  first  locates  the 
inboard  and  outboard  ends  of  the  slats  and  defines  the  slat  length.  The  inboard 
location  is  set  at  1.  5 ft.  outboard  of  the  side  of  the  body.  The  outboard  location 
includes  3.0  ft.  of  clearance  for  each  wing  mounted  engine  pylon.  The  second  opera- 
tion determines  the  individual  slat  lengths,  chords,  and  inboard  and  outboard  stations 
for  two  and  four  engine  aircraft.  The  slat  analysi  5 for  a two-engine  configuration 
provides  three  options  for  slat  segment  location-  1)  inboard  only,  2)  outboard  only,  3) 
outboard  cnly,  4)  inboard  and  center,  5)  center  and  outboard,  and  6)  inboard,  center, 
and  outboard.  The  specific  slat  chord  lengths  are  computed  as  a function  of  the  slat 
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chord  to  wing  chord  ratio.  However,  if  the  ratio  is  not  input  a value  of  0.  0735  is 
used.  This  is  an  average  value  for  typical  transport  aircraft. 

The  fixed  leading  edge  geometry  routines  provide  planform  dimensions  and  locations 
for  the  wing,  horizontal  stabilizer,  and  vertical  stabilizer  leading  edges.  The 
horizontal  and  vertical  stabilizer  leading  edges  start  at  the  body  and  end  at  the  tip. 
The  leading  edge  chord  is  input  as  the  total  distance  forward  of  the  front  spar.  The 
wing  has  two  types  of  fixed  leading  edges;  under-slat  and  be*vveen-slat.  The  leading 
edges  extend  from  the  side  of  the  body  to  the  tip,  the  appropriate  type  being  used  as 
a function  of  the  slat  locations.  The  between-slat  type  extends  the  full  distance 
forward  of  the  front  spar  and  the  under-slat  type  assumes  a chord  equal  to  8^  of  the 
wing  chord.  Fixed  leading  edge  geometry  output  consists  of  the  lengths  and  chords  of 
each  type  of  edge. 

The  spoiler  geometry  routine  provides  spoiler  p!  inform  dimensions  and  locations 
from  input  data.  If  the  spoiler  area  is  input  the  . iler  will  be  sized  to  the  area 
output  from  the  aircraft  sizing  routine.  If  the  area  is  not  input  the  user  must  provide 
the  inboard  and  outboard  edge  locations  as  well  as  the  spoiler  chord  to  wing  chord 
ratio.  If  the  spoiler  chord  to  wing  chord  ratio  is  not  input  it  is  assumed  to  be  0.  15. 
The  spoiler  inboard  edge  is  assumed  to  be  at  the  side  of  the  body  and  the  outboard 
edge  is  computed.  The  outboard  edge  is  truncated  at  the  wing  tip  or  at  the  edge  of 
the  aileron.  Spoiler  geometry  output  consissts  of  inboard  and  outboard  chords,  span 
stations  of  spoiler  inboard  and  outboard  edges,  and  the  spoiler  length. 

The  fixed  trailing  edge  geometry  routine  assumes  a total  length  from  the  body  to  the 
tip  tor  wings,  horizontal  stabilizers,  and  vertical  stabilizers.  The  fixed  trailing 
edge  chord  is  computed  as  a function  of  the  total  trailing  edge  length  and  the  sur- 
faces involved.  The  lower  surface  chord  is  computed  as  G.  8%  of  the  trailing  edge 
length  if  there  are  flaps  and  I0cl  if  there  are  ailerons,  ruddei  c or  elevators.  The 
upper  surface  chord  is  computed  as  29.  Gfi  of  the  trailing  edge  length  for  flaps  only. 

It  is  set  equal  to  the  spoiler  chord  if  there  are  flaps  and  spoilers,  and  equal  to  10ca 
of  the  trailing  edge  length  for  ailerons,  rudder,  or  elevators.  If  there  are  no 
control  surfaces  the  fixed  trailing  edge  extends  from  the  rear  spar  to  the  aft  edge  of 
the  wing,  horizontal  stabilizer,  or  vertical  stabilizer. 

The  spoiler  analysis  produces  structural  member  thicknesses  and  desired  rivet  pat- 
terns. The  planform  geomcti’y  is  obtained  from  the  spoiler  geometry  output. 

Member  thicknesses  are  computed  ;ind  adjusted  to  standard  gages.  Cross-sectional 
geometry  is  shown  in  Figure  16.  The  front  spar  is  a bent-up  sheet  metal  zee,  the 
two  ribs  (at  each  support)  are  bent-up  sheet,  and  the  skins  are  sheet  metal  over  a 
full  depth  honeycomb  core. 
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The  spoiler  analysis  accounts  for  external  and  internal  loads.  The  external  loads 
for  transport  aircraft  are  normally  those  loads  which  the  spoiler  actuator  produces. 
In  this  analysis  the  spoiler  external  load  is  assumed  to  be  600  lb  - in.  of  hinge 
torque  per  running  inch,  limit.  This  is  comparable  to  the  Model  990  loading  con- 
dition. The  internal  load  analysis  subdivides  the  total  spoiler  area  into  the  smallest 
number  of  segments  (individual  surfaces)  where  all  segments  are  equal  in  length  and 
not  longer  than  60  inches.  The  segments  are  supported  at  each  end  and  all  torque  is 
taken  by  the  inboard  support.  The  spoiler  is  analyzed  as  a simple  beam.  The  point 
of  maximum  bending  moment  is  determined,  and  the  bending  moment  and  spar  depth 
computed.  All  spoiler  bending  moment  is  taken  by  the  spar  and  effective  skin.  The 
bending  section  (Figure  16)  is  assumed  symmetrical,  and  the  tension  and  compres- 
sion stresses  are  equal  to: 

M(d/2) 


where 


bending  stress 
contour  depth  at  spar 
bending  moment 
section  moment  of  inertia 


The  compression  buckling  allowable  is 


0.56  F 


cy  \ A 


(Reference  8,  Eq. 
C7.  4) 


whe  re 


F = compression  buckling  allowable 
c s 

F = compressive  yield  allowable 

t = material  thickness 

A = cap  area  (=  1 . 73t) 

E = material  elastic  modulus 


The  spar  cap  sheet  thickness  is  sized  so  that  the  stress  level  is  equal  to  or  less  than 
the  larger  of  the  compression  buckling  allowable  or  S 0 Vo  of  the  ultimate  tensile  allow- 
able. 


The  ir.hoard  rib  is  analyzed  for  bending  at  the  front  spar.  Since  all  torque  is  taken  at 
this  rib,  the  bending  moment  is  equal  to  the  total  spoiler  torque  about  the  spar.  The 
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section  (Figure  16)  is  symmetrical  and  the  tension,  compression,  and  compression 
buckling  stresses  are  computed  the  same  as  shown  for  the  spar. 

The  skin  thickness  is  based  on  skin  shear  flow  at  the  inboard  hinge  where  all  spoiler 
torque  is  reacted  about  the  spar.  Since  the  skin  is  supported  by  the  honeycomb  core, 
the  shear  allowable  is  based  on  the  ultimate  shear  stress  times  a rivet  factor  of  0.  8. 

Appropriate  material  properties  are  selected  for  each  part  analyzed.  The  analysis 
determines  the  material  thicknesses  as  a minimum  required  thicknesses  and  then 
rounds  the  value  of  the  next  larger  standard  gage.  A minimum  gage  of  0.  020  in.  and 
a maximum  gage  of  0.250  in.  are  set  as  constraints.  The  standard  sheet  gages  used 
are  summarized  in  Table  1. 


Table  1.  Standard  Sneet  Gages 


in. 

in. 

0.020t 

0.071 

0.  025 

0.080 

0.032 

0.090 

0.036 

0. 100 

0.040 

0.125 

0.045 

0.160 

0.  050 

0.190 

0.063 

0.250# 

t Minimum 

# Maximum 

The  number  of  rivet  holes  (repre- 
senting the  actual  number  of  rivets 
needed)  and  the  hole  sizes  are  out- 
put. The  quantity  and  size  of  the 
rivets  is  based  on  a T/2A  shear 
flow  analysis  at  the  inboard  rib. 

The  rivets  are  sized  based  on  the 
protruding  head  shear  allowables 
at  a spacing  of  four  times  the 
shank  diameter.  The  number  of 
holes  is  equal  to  the  number  of 
rivets.  That  is,  the  holes  are 
counted  for  only  one  member.  When 
two  rows  of  rivets  are  required, 
an  additional  amount  of  spar  or  rib 
cap  width  is  output,  but  the  addition- 
al area  is  not  used  to  resize  the  cap. 


The  foreflap  analysis  produces  the  structural  member  dimensions  and  desired  rivet 
patterns.  The  planform  geometry  is  obtained  from  the  foreflap  geometry  output. 
Member  thicknesses  are  computed  and  then  adjusted  to  standard  gages.  A typical 
foreflap  cross  section  is  shown  in  Figure  17.  The  front  spar  is  bent -up  sheet  metal 
channel  and  is  sized  by  a loads  analysis.  The  leading  edge  skin  and  rib  thicknesses 
are  fixed  at  0.  050  in.  The  honeycomb  box  factor  is  set  at  1 and  assumes  an  allow- 
able shear  stress  of  160  psi.  The  box  skin  thickness  is  assumed  to  be  0.020  in. 

The  foreflap  spar  analysis  accounts  for  external  and  internal  loads.  The  external 
applied  loads  are  derived  from  the  general  formula: 


W = S C — 
n \ 295 


Figure  17.  Foreflap  Geometry 


total  surface  load 


total  surface  area 


normal  lift  coefficient 


design  speed 


The  average  pressure,  ultimate,  is  applied  to  the  foreflap  uniformly  and  is  computed 
from  the  transposed  form: 


ave  S \ 295 


where 


P = average  ultimate  surface  pressure  and  for  the  foreflap 
ave 

C = 4.0 
n 


The  internal  load  analysis  subdivides  the  total  surface  length  into  a number  of  equal 
length  segments  (individual  surfaces)  each  with  a length  equal  to  or  less  than  180  in. 
If  the  individual  segment  length  turns  out  to  be  greater  than  140  in.  , three  hinge 
supports  are  assumed.  One  is  in  the  center  and  two  are  located  15%  of  the  surface 
length  from  each  end.  If  the  individual  surface  length  is  less  than  or  equal  to  140  in. 
two  hinge  supports  are  assumed,  each  28%  of  the  surface  length  from  each  end. 


The  vertical  shear,  bending  moment,  and  torque  about  the  front  spar  are  calculated 
at  each  hinge.  The  torque  is  calculated  at  each  end  of  the  surface  segment  and  is 
assumed  to  vary  linearly  between  the  ends.  T^e  torque  is  reacted  at  each  hinge 
using  the  same  formulae  used  to  calculate  shear  reactions.  The  foreflap  bending  is 
assumed  to  be  taken  by  the  spar  and  associated  skin  as  shown  in  Figure  17.  The 
bending  stress  can  be  computed  from  Equation  1,  and  the  compression  buckling 
allowable  stress  can  lie  computed  from  Equation  2.  Spar  thickness  is  sized  to  be  the 
minimum  necessary  so  that  the  stress  level  is  equal  to  or  less  than  the  larger  of  the 
compression  budding  allowable  or  80%  of  the  ultimate  tensile  allowable. 


All  rivet  patterns  are  assumed  to  be  comprised  of  a single  row  of  0.  25  in.  diameter 
rivets  spaced  at  two  diameters.  The  output,  number  of  holes  is  equal  to  the  number 
of  rivets.  However,  each  rivet  is  accounted  for  in  only  one  part  of  the  joint. 
Adjustment  of  material  thicknesses  to  a standard  gage  is  accomplished  in  the  same 
manner  as  discussed  for  the  spoiler. 
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The  analysis  of  the  flaps,  ailerons,  slats,  rudder,  and  elevators  produces  the 
structural  member  dimensions  and  the  desired  rivet  patterns.  The  planform  geo- 
metry is  obtained  from  the  specific  control  surface  geometry  output,  and  the  mem!  er 
thicknesses  are  computed  and  then  adjusted  to  standard  gages.  The  control  surfaces 
are  assumed  to  have  the  geometry  shown  in  Figure  18.  The  front  spar  has  extended 
caps  and  a sheet  metal  web,  and  the  rear  spar  is  a bent-up  sheet.  Both  the  leading 
edge  skin  and  the  main  box  skin  are  sheet  metal.  The  trailing  edge  consists  of  a 
full -depth  honeycomb  core  with  a single  piece  of  sheet  metal  forming  both  upper  and 
lower  skins.  The  airload  ribs  and  the  leading  edge  ribs  are  bent-up  metal.  There  is 
a leading  edge  rib  at  each  airload  rib  span  station.  The  hinge  ribs  consist  of  ex- 
truded spar  caps  and  a sheet  metal  web  with  bent-up  flanges  to  pick  up  front  and 
rear  spars. 

Appropriate  material  properties  are  selected  for  the  analysis  of  each  part.  Thick- 
nesses are  fixed  for  the  leading  edge  skin  and  ribs,  airload  ribs,  rear  spar,  and 
trailing  edge  skin  as  follows: 


Leading  edge  skin 
Leading  edge  ribs 
Airload  ribs 
Rear  spar 
Trailing  edge  skin 


Thickness 

Same  as  box  skin 
Same  as  airload  ribs 
One  gage  heavier  than  skin 
One  gage  heavier  than  skin 
Minimum  gage 


The  analysis  for  the  remaining  parts  determines  the  material  thicknesses  in  terms  of 
a minimum  required  thickness  and  then  rounds  the  value  to  the  next  larger  standard 

gage.  Standard  sheet  gages  are  summarized 
Table  2.  Standard  Extrusion  Gages  in  Table  1,  and  standard  gages  for  extru- 

sions  in  Table  2. 


~ 7T7T The  parts  sized  by  a loads  analysis  include 

q ^ 0 156  the  °asic  s^lns>  sPar  webs,  spar  caps, 

, ' hinge  rib  caps,  hinge  rib  webs,  and  the 

0.078  0.188  , b,.  , ’ , „ ’ . 

q ^ q trailing  edge  honeycomb.  The  analysis 

’ ’ accounts  for  both  the  internal  and  external 

4,  . ,,  . loading  conditions.  The  applied  external 

T Minimum  # Maximum  , , , 

loads  are  normal  (to  the  surface)  loads 

only.  For  the  wing  surfaces  (flaps,  ailerons,  and  slats)  these  normal  loads  are  de- 
rived from  the  general  formulae  Equations  3 and  4. 


For  flaps, 


V = 1.75  V (Ref.  MIL-8860.  Para.  6.2.3.  9),  where  V = stall  speed 
s s 


FRONT 
SPAR 
0.  15  c 


For  slats, 


V 

C 


1.75  V 


3.0 


For  ailerons,  rudders,  and  elevators,  V is  derived  from 

2 


V 

N W = C T S — 

z L , Wing  295 
Max 


(MIL-8860,  Para.  3. 2. 2.  2); 


or  transposing: 


V = V = — 
a C 


295  N W 
z 


whe  re 


N 


L Wing 
Max 


maximum  normal  load  factor 


W 

C. 


Max 


^Wing 


V 


For  ailerons, 


C =1.6 
n 


aircraft  gross  weight 
maximum  lift  coefficient 

wdng  area 

aileron  design  speed 


For  rudders  and  elevators , 


C =1.3 

n 

The  average  pressure,  P , is  applied  to  the  control  surface  as  a chordwise  triangu- 
lar distribution  with  the  center  of  pressure  at  the  33^  chord  aft  of  the  leading  edge. 

If  the  design  speed  is  equal  to  or  greater  than  Mach  1,  the  center  of  pressure  for  the 
aileron,  rudder,  or  elevator  is  assumed  to  be  at  the  47f,"o  surface  choi’d.  Spanwise 
running  surface  loads  are  there  fore  proportional  to  surface  chord. 

The  internal  load  analysis  subdivides  the  total  surface  length  into  a number  of  equal 
length  segments  (individual  surfaces)  each  with  a length  equal  to  or  less  than  180  in. 
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If  the  individual  -segment  length  is  140  in.  or  less,  two  hinge  supports  are  assumed, 
located  28%  of  the  total  length  from  each  end.  If  the  segments  are  greater  than  140 
in.  , three  hinge  supports  are  assumed.  One  is  located  in  the  center  and  two  are 
located  15%  of  the  total  length  from  each  end.  The  vertical  shear,  bending  moment, 
and  torque  about  the  front  spar  are  claculated  at  each  hinge.  Torque  is  calculated  at 
each  end  of  the  surface  segment  and  Is  assumed  to  vary  linearly  between  the  ends. 

For  flaps  and  slats,  torque  is  reacted  at  each  hinge  using  the  same  formulae  used 
to  calculate  shear  reactions.  For  ailerons,  rudders,  and  elevators  all  torque  is 
reacted  at  the  inboard  (or  lower)  hinge. 

The  skin  thickness  is  computed  based  on  skin  shear  flow,  and  the  allowable  stresses 
are  fixed  as  a function  of  rib  spacing.  Since  the  hinge  rib  number  and  locations  are 
fixed,  rib  spacing  is  determined  for  each  bay  between  hinge  ribs  by  equally  spacing 
airload  ribs.  For  a given  skin  thickness,  rib  spacing  can  be  determined  from 
Figure  19.  This  curve  is  a typical  design  curve  for  sonic  fatigue  requirements  as 
discussed  in  Lockheed  Design  Handbook,  pages  10.801  - 10.809)  and  Lockheed  Report 
LG1538-3-1,  Figure  5.  3. 


figure  19.  Sonic  Fatigue  Curve. 


t 

(• 
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For  practical  considerations,  a minimum  rib  spacing  of  3.0  in.  is  used. 


An  analysis  is  made  of  the  inboard  panel  of  the  bay  with  the  maximum  rib  spacing 
assuming  maximum  skin  shear  flow  exists  there.  Allowables  are  determined  for  an 
incomplete-diagonal -tension  panel  utilizing  I^ACA  TN  2661.  The  critical  buckling 


stress  is  computed  from  F = K E (t/d)“  where  K is  from  Figure  12  of 

ss  ss 

CR 


Reference  2.  The  diagonal  tension  factor,  K,  is  derived  from  Equation  27  of 
Reference  2.  Then  the  allowable  shear  stress  can  be  determined  as  a function  of  K 
utilizing  the  40-degree  curve  of  Figure  19  (a)  of  Reference  2.  'Hie  skin  is  sized  so 
that  the  maximum  shear  stress  does  not  exceed  the  allowable,  and  so  that  the  ratio 
of  the  maximum  to  the  critical  shear  stress  does  not  exceed  5. 


The  spar  web  thickness  is  determined  using  the  maximum  spar  shear  flow.  The 
analysis  is  made  using  either  the  panel  at  the  in  bo  a I'd  end  of  the  surface  segment  or 
the  panel  just  outboard  of  the  inboard  hinge,  which  ever  has  the  greatest  ratio  of 
spar  height  to  rib  spacing.  An  incomplete  diagonal -tension  analysis  is  made  like 
that  made  for  the  skin. 


All  flap  bending  moment  is  taken  by  the  front  spar  caps  and  associated  skin  and  spar 
web.  The  critical  bending  location  is  at  the  hinge  where  the  ratio  of  bending  moment- 
to-spar  depth  is  largest.  The  effective  spar  section  is  as  shown  in  Figure  IB. 


This  bending  section  is  symmetric;  therefore,  tension  and  eompresrion  stresses  are 
equal  and  may  be  computed  from  Equation  1. 


d = contour  depth  at  spar 


The  compression  buckling  allowable, 


Ts  J(Tj(r)'  ’ 


(Reference  S , Equation  C7.5i 


whe  re 


F compression  buckling  allowable 
cs 


F = compressive  vield  allowable 

cy 


t = material  thickness 


!•  ’ 


k 


■j  ),iui  I'fir.liiiwiuyiipy^WllpiljpjWgawwiWf^-W 


A = cap  area  (=  1.  46  t-  + 0.  82t) 
E = material  elastic  modulus 


The  spar  cap  is  assumed  to  be  an  extrusion  with  a constant  section  thickness  sized 
so  that  the  stress  level  is  equal  to  or  less  than  the  larger  of  the  compression  buckling 
allowable  or  80%  of  the  ultimate  tensile  allowable. 

For  all  surface  types,  hinge  ribs  are  assumed  to  have  the  same  part  thickness  as  the 
inboard  hinge.  The  rib  cap  is  sized  by  the  rib  bending  moment  at  the  front  spar 
which  is  equal  to  the  surface  torque  (about  front  spar)  at  the  inboard  hinge.  The 
generalized  effective  rib  section  is  considered  to  be  the  same  as  the  spar  section. 

The  compression  buckling  allowable  stress  equation  is  the  same  as  that  used  for  the 
spar.  The  rib  cap  is  assumed  to  be  an  extrusion,  and  the  constant  section  thickness 
is  sized  in  the  same  manner  as  the  spar  cap.  The  web  thickness  is  sized  to  be 


adequate  for  the  inboard  hinge  rib  shear  flow, 


where 


Q = inboard  hinge  shear  flow 


T = torque  reacted  by  the  inboard  hinge 
A = inter-spar  box  area  at  the  inboard  hinge 

The  shear  buckling  stress  is  calculated  for  a web  panel  at  the  front  spar  assuming  a 
panel  aspect  ratio  of  2. 


SCR 


shear  buckling  stress 
material  elastic  modulus 
material  thi  cknes  s 
front  spar  height  at  rib 


The  web  thickness  is  sized  so  that  the  shear  stress  level  is  equal  to  or  less  than  the 
larger  of  the  shear  buckling  stress  or  80%  of  the  ultimate  shear  allowable. 
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The  assumed  honeycomb  type  and  size  has  a shear  allowable  ol  160  lb/in  . A factor 
is  developed  that  indicates  how  much  heavier  than  the  basic  core  the  actual  core 
must  be. 
load. 


The  factor,  K 


core 


is  based  on  the  core  shear  due  to  trailing  edge  air- 


1 


(0. 2 P ) (0. 2 chord) 
max 


where 


- core  shear 


= maximum  airload 


chord  = chord  length 


d = contour  height  at  rear  spar 


K = F / 1 60 
core  s 


Rivet  sizes  and  numbers  are  calculated  using  the  shei  \ flows  that  sized  the  skin, 
spar  web,  and  hinge  rib  web.  Rivet  shear  values  are  used  as  the  allowables  and  a 
rivet  spacing  of  four  diameters  is  assumed. 


(Ib/in) 


0 to  776 
777  to  954 
955  to  1573 
1574  to  1957 
1908  and  above 


No.  of 
Rows 


Spacing 


Rivet 


0.  50 
0.  625 
0.75 
0.625 
0.75 


In  the  output  the  number  of  holes  is  equal  to  the  number  of  rivets;  each  rivet  hole  is 
accounted  for  in  only  one  part  of  the  joint.  When  two  rows  of  rivets  are  required,  an 
additional  spar  or  rib  cap  width  is  output.  This  additional  area  is  not  used  to  re- 
size the  cap. 


3.2  TIP,  LEADING  AND  TRAILING  EDGE  PART  DEFINITION 

The  tip,  leading  edge,  and  trailing  edge  part  definition  routines  define  the  detail 
parts  making  up  the  fixed  leading  edge,  fixed  trailing  edge,  slats,  flaps,  foreflaps, 
control  surfaces  (spoilers,  ailerons,  rudder,  and  eleva;ors) , and  tips.  The  data 
that  is  generated  includes  the  number  of  parts,  part  dimensions,  weight,  and  cost 
parameters.  The  parts  definition  derives  its  input  from  previous  geometry  and 
analysis  subroutines. 


The  fixed  leading  edge  segments,  as  defined  by  the  geometry  subprogram,  are 
divided  into  a number  of  60-in.  sections  with  one  shorter  section.  If  the  segment 
is  60  in.  or  less,  only  one  section  is  assumed.  The  under-slat  leading  edge  is  made 
of  two  skins  splice  J at  the  nose  with  an  extruded  angle  (chafing  strip).  The  between- 
slat  leading  edge  has  a one  piece  skin;  the  skin  perimeter  is  assumed  equal  to  2.  5 
times  the  fixed  leading  edge  chord  The  upper  skin  of  the  under-slat  segment 
utilizes  a factor  of  1 . 5 and  the  lower  skin  a factor  of  1.0.  The  skin  thickness  is  set 
at  0.040  in.  with  the  chafing  strips  and  edge  member  thicknesses  set  at  0.  060  in. 

The  ribs  are  spaced  at  10-in.  increments,  and  the  rib  height  is  assumed  to  be  0.  85 
times  the  rib  chord  length,  The  ribs  are  made  of  bent-up  0.040-in  sheet  with 
lightening  holes.  The  rib-to-skin  fasteners  are  5/32-in.  diameter  rivets  spaced  at 
0.75-in  intervals.  The  chafing  strip  rivers  are  5/32  in.  in  diameter  spaced  at 
0.625-in.  intervals,  and  the  edge  member-to-skin  rivets  are  3/i6  in.  in  diameter 
spaced  at  1.5  in. 

llie  fixed  trailing  edges  for  the  wings,  horizontal  stabilizer,  and  vertical  stabilizer, 
illustrated  in  Figure  20,  are  assumed  to  be  comprised  of  flat  sheet  skins  and  bent- 
up  sheet  ribs.  All  skins  are  0.037-in.  thick  and,  like  the  fixed  leading  edge,  are 
defined  in  terms  of  60-in.  segments.  The  ribs  are  spaced  at  10-in.  increments  and 
are  constructed  of  bent -up  0.  040-in.  sheet  with  a 0.  73-in.  flange  on  each  edge. 
Lightening  holes  are  spaced  at  1. 5-in.  intervals  and  have  a diameter  of  0.  375  times 
the  local  chord.  The  skins  attach  along  the  forward  edge  and  along  each  rib  with 
5/32-in.  diameter  rivets  spaced  at  four  diameters. 

The  spoiler,  illustrated  in  Figure  16,  is  assumed  to  be  comprised  of  a spar,  skins, 
honeycomb  core,  and  a wedge  shaped  skin  elosui'e.  The  parts  definition  process 
defines  the  dimensions,  and  the  rivet  sizes  and  quantities  based  on  the  spoiler 
stress  analysis.  The  material  weight  assumes  1.0  in.  added  to  the  length  and  width 
dimensions  of  the  sheet  Hat  pattern,  and  to  all  dimensions  of  tl^e  full-depth  honey- 
comb core.  The  material  weight  of  the  core  includes  0. 1 lb/ft  for  adhesive. 

The  parts  definition  for  the  foreflap  (Figure  17)  derives  the  dimensions,  end  the  rivet 
sizes  and  quantities  from  the  foreflap  stress  analysis.  The  upper,  lower,  and 
leading  edge  skins  have  material  weights  calculated  assuming  1.0  in.  of  additional 
material  on  all  sides.  The  leading  edge  skin  width,  or  cross-section  periphery,  is 
set  equal  to  2.  64  times  leading  edge  chord.  Foreflap  cross-sectional  area  aft  of  the 
spar  is  calculated  as 

Area  - (spar  height)  - (chord  length  aft  of  spar)  ' 0.698 


This  formula  provides  the  basis  for  computing  the  honeycomb  core  and  closing  rib 
weights.  Material  weight  for  the  core  is  based  on  maximum  dimensions  plus  1.  0 in. 
Closing  rib  material  weight  is  based  on  Hat  pattern  dimensions  plus  1. 0 in.  on  each 
side. 


C = LESSOR  OF  B OR  D/2 
Cl  = LESSOR  OF  B OR  D/2 


LOCATION 

B 

bl 

AT  FLAP  (NO  SPOILER) 

0.  296A 

0.  068A 

AT  FLAP  (INBD./OUTBD. 

SPOILER 

0.  068A 

OF  SPOILER) 

CHORD 

AT  AILERON 

0.  10A 

0.  10A 

AT  RUDDER/ELEVATOR 

0.  10A 

0.  10A 

INBOARD  OF  AILERON  (NO  FLAP) 


Figure  20  Fixed  Trailing  Edge. 


The  parts  definition  process  tor  the  flaps,  ailerons,  rudders,  elevators,  and  slats 
(Figure  18)  derives  the  dimensions,  and  rivet  sizes  and  quantities  from  the  control 
surface  stress  analysis.  The  surface  skins  are  assumed  to  be  made  in  three  pieces. 
The  inboard  and  outboard  skins  are  assumed  to  have  a length  equal  to  28%  of  the  sur- 
face length  and  the  center  44%  of  the  surface  length.  The  leading  edge  skin  width 
(periphery  of  leading  edge  cross  section)  is  calculated  from  the  following: 


inboard  skin  .,idth.  INSWI  - K j ^DCSWI)  - . 28  (DCgWI.- DCSWO) 
Center  skin  width,  DN?WC  = K 


(.  15) 


2 (DCSWI)  - (DCSWI  - DCS  WO) 
2 


(.15) 


DCSWO  + DCSWI  - .72  (DCSWI  - DCSWO) 


(.15) 


Outboard  skin  width,  DNSWO  = K 
where 

K =2.98  for  slats 

= 2. 57  for  other  surfaces 
DCSWI  = inboard  chord  length  of  surface 
DCSWO  = outboard  chord  length  oi  surface 


Computation  of  the  front  spar  and  hinge  rib  cap  material  weight  assumes  an  additional 
2.0  in.  on  the  extrusion  length.  Rear  spar  material  weight  assumes  an  additional 
0.  5 in.  cn  all  sides  of  the  flap  pattern  dimensions.  Material  weight  for  the  skins  is 
computed  as  the  actual  weight  plus  0. 5 in.  of  additional  material  on  all  edges.  Of 
the  total  skin  rivets  32%  are  assumed  to  be  in  each  of  the  inboard  and  outboard  skins, 
and  36%  in  the  center  skin. 

Airload  ribs  are  bent-up  sheet  metal  and  material  weight  is  based  on  the  flat  pattern 
dimensions  plus  an  additional  1.  0 in.  in  both  length  and  width.  Theoretical  and 
actual  rib  weights  assume  lightening  holes  with  diameter  equal  to  75%  of  average  rib 
height  spaced  at  1-1/2  diameters. 

The  nose  ribs  are  assumed  to  be  parabolic.  Material  weight  is  based  on  1.0  in. 
added  to  the  length  and  width  of  the  flat  pattern  dimensions.  Each  rib  contains  one 
lightening  hole  with  a diameter  equal  to  7 5%  of  the  smaller  rib  chord  length  or  84.  5 % 
of  rib  height.  The  hinge  rib  webs  are  a solid  web  with  no  lightening  holes.  Material 
weight  is  calculated  assuming  0.  5 in.  of  additional  material  on  all  edges.  The 
honeycomb  trailing  edge  wedge  theoretical  weight  is  computed  as  the  theoretical 
weight  times  the  honeycomb  core  factor  from  the  stress  analysis  routines.  Material 
weight  is  computed  assuming  a honeycomb  block  with  dimensions  equalling  the 
largest  web  dimensions  plus  1.0  in.  and  adhesive  weight. 
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The  parts  definition  process  for  the  tip  assumes  the  geometry  and  part  dimensions 
shown  in  Figure  21.  Actual  weight  for  the  skin  is  computed  from: 


WT  = 30  (0.032)  (TIP  CHORD;  (DENSITY) 


The  material  weight  for  all  sheet  metal  pails  assumes  an  additional  1.0  in.  of 
material  on  both  the  length  and  the  width.  All  attachments  assume  a single  row  of 
3/l6-in.  diameter  rivets  spaced  at  four  diameters 


3 . 3 COM  PU  T E R PROG  RAM 


The  computer  program  is  a series  of  modularized  subroutines,  developed  for  the 
CDC  6400  computer  and  now  operational  on  the  6600,  to  define  geometry,  perform 
structural  analysis,  and  develop  part  definitions  for  aerodynamic  surface  leading  edge, 
trailing  edge  and  tips.  A listing  of  these  subroutines  is  separately  available. 


The  subroutines  have  been  developed  in  such  a manner  to  allow  eas}  integration  of 
any  or  all  subroutines  within  existing  or  future  developmait  of  aircraft  sizing  and 
cost  synthesis  programs.  To  provide  this  end  the  following  ground  rules  were 
employed: 


Total  modularity. 

Control  subroutine  with  input. 

Design  for  ease  of  update,  modification  and  integration. 
Level  of  analysis  is  function  of  input. 


3.3. 1 PROGRAM  DESCRIPTION.  The  program  as  described  here  is  a series  of 
individual  modular  programs  that  operate  as  individual  entities.  The  generalized 
flow  of  information  from  one  subroutine  to  another  is  as  shown  in  Figure  ’5.  However, 
each  subroutine  provides  only  one  part  of  the  total  information  required  for  the  over- 
all system. 


The  computer  program  was  designed  in  a modular  fashion  as  shown  in  tire  general 
overlay  structure  of  Figure  22.  The  driver  has  supporting  subroutines  resident  and  tire 
the  primary  and  secondary  overlays  are  divided  bv  function.  The  functions  include: 
control  surfaces,  tips,  fixed  leading  edge,  spoilers,  fixed  trailing  edge,  and  geo- 
metry. This  overlay  provides  an  operable  computer  program  for  easy  integration 
into  a total  synthesis  operation. 


The  subroutine  core  requirements,  based  on  the  overlay  structure  shown  in  figure 
21  is  outlined  in  luble  3.  These  core  requirements  are  the  maximum  for  each 
operation. 


Figure  22.  Typical  Overlay  Structure. 
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Table  3.  Subroutine  Core  Requirements. 


Overlay  Designation 


Core  Required  (Octal) 


DRV  + 1,0  + 1.1 
DRV  + 1,0  + 1,2 
DRV  +1,0  +1,3 
DRV  + 1,0  + 1,4 


42500 

35000 

33500 

32500 


DRV  +2,0 


31100 


DRV  +3,0 


34000 


DRV  + 4,0 


36000 


DRV  + 5,0 


35300 


DRV  +6,0 


35000 


3.3.2  DESCRIPTION  OF  SUBROUTINES.  The  subroutines  that  were  developed 
within  this  study  effort  are  described  as  follows: 


DRV  is  the  driver  subroutine  that  acts  as  the  main  control  routine  for  the  total  pro- 
gram. It  calls  the  subroutine  in  the  proper  logical  sequence  as  required  to  perform 
the  task. 


INITIAL  initializes  all  variables  to  their  required  predetermined  value  and  sets  all 
input  variables  to  zero. 


THICK  defines  control  surface  cross-section  depth  as  a function  of  type,  span,  and 
chord  location. 


CHANGE  selects  the  proper  material  cone  for  each  part  as  a function  of  input  flags 


DISC  Stores  part  sizes  and  cost  parameters. 


PROPMT  stores  the  mechanical  properties  of  various  materials  for  use  in  the  stress 
analysis  and  part  definition  routines. 


SMRQOT  calculates  the  root  of  a given  quadratic  eqution  as  required  by  the  different 
geometry  routines. 


DENS  calls  the  change  routine  and  provides  the  appropriate  material  densities  to  the 
part  definition  routines. 
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CALLSF  subdivides  the  control  surfaces  into  equal  length  segment,  no  longer  than 
180  inches.  It  also  provides  the  communication  link  with  the  appropriate  analysis 
and  parts  definition  routines  for  each  segment. 


A SURF  provides  the  structural  analysis  for  the  flaps,  ailerons,  slats,  rudder,  and 
elevators.  This  routine  provides  part  dimensions,  actual  ;uid  theoretical  thicknesses, 
and  rivet  sizes  and  quantities. 


CSPART  ut'.izes  data  from  ASURF  to  produce  the  number,  sizes,  weights,  and  cost 
parameters  of  all  parts  for  flaps,  ailerons,  rudders,  elevators,  and  slats. 


FORFLP  computes  spar  thicknesses  and  rivet  patterns  for  the  loreflap.  The 
structural  analysis  is  provided  by  ASURF  through  the  flap  routine. 


FFLPPT  utilizes  data  from  FORFLP  to  produce  the  number,  sizes,  weights  and 
cost  parameters  for  all  parts  of  the  foreflap. 


T1PART  produces  the  number,  sizes,  weight,  and  cost  parameters  for  all  surface 
tip  parts. 


FLPART  utilizes  data  from  FLEGOM  to  select  planform  dimensions  for  each 
leading  eage  segment  and  calls  the  part  definition  routine  {FXLEPTj  for  that  seg- 
ment. 


PXLEPT  produces  the  number,  sizes,  weight,  mid  cost  parameters  for  all  leading 
edge  parts. 


CALLSP  provides  the  spoiler  into  equal  segments  of  GO-inch  lengths  maximum.  It 
also  serves  as  the  calling  routine  for  spoiler  analysis  mid  spoiler  and  spoiler  part 
definition  routines. 


SPOIL  provides  the  structural  analysis  for  the  spoilers.  This  routine  provides  part 
dimensions,  actual  and  theoretical  thickness,  :v>d  rivet  sizes  mid  quantities. 


SPPART  utilizes  data  from  the  spoil  routine  to  produce  the  number,  sizes,  weights, 
and  cost  parameters  of  all  parts  for  the  spoilers. 


FEPART  provides  quantity,  sizes,  weights,  and  cost  parameters  for  fixed  trailing 
edge  parts  on  wings,  horizontal  stabilizers,  and  vertical  fins. 


SLATLO  locates  inboard  and  outboard  ends  of  slats  for  use  in  the  SLGEOM  routine. 


SLGEQM  defines  slat  planform  dimensions  and  locations  from  input  data. 


EVGEOM  defines  elevator  planform  dimensions  and  locations  from  input  data. 


RUDGOM  defines  rudder  planform  dimensions  mid  locations  from  input  data. 


AILGOM  defines  aileron  planform  dimensions  mid  locations  from  input  data. 


FLPGOM  defines  foreflap  and  flap  planform  dimensions  and  locations  from  input  data. 


SPLGOM  defines  spoiler  planform  dimensions  and  locations  from  input  data. 


FLEGOM  defines  fixed  leading  edge  planform  dimensions  and  locations  for  wings 
horizontal  stabilizers,  and  vertical  fins. 
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3.3.3  USER'S  GUIDE.  This  computer  program , LTTWPD,  was  developed  on  a 
CDC  6400  computer.  The  Convair  Aerospace  Computer  Program  Library  No.  is 
P5749  and  a listing  of  the  program  is  separately  available.  This  program  was 
written  in  Fortran  IV  language  in  an  update  file  format  It  is  stored  on  magnetic 
tape  in  the  Convair  Tape  Library,  Tape  No.  14289. 

A namelist  type  input  was  used  with  the  variable  name,  value,  and  definition  of 
variable  name  on  one  card.  A complete  list  of  input  variables  is  presented  in 
Table  4.  The  deck  setup  for  executing  the  program  is  illustrated  in  Figure  23. 

The  input  variables  required  for  each  structural  type  are  presented  in  sample 
case  listings  (Tables  5 through  13).  All  input  variables  are  initialized  to  zero  for 
each  case.  Therefore,  for  multiple  case  runs  you  must  input  variables  for  each 
case. 

The  output  consists  of  input,  detail  weight  and  part  definition  data,  and  total 
weight.  Sample  cases  are  presented  in  Table  14. 


Table  4.  List  of  Input. 


PSINPNL 

c 

c 

c 


CLESI  - 

0.0 

CLESO  ■ 

0.0 

CLMAX  ■ 

0.0 

CSCDWC- 

0.0 

CW5PI N» 

0.0 

CW5POT- 

0.0 

0 ■ 

0.0 

DBW  » 

0.0 

DCSL  - 

0.0 

OCSWI  - 

0.0 

DCSWO  ■ 

0.0 

DCSYI  - 

0.0 

ocrsi  * 

0.0 

DCYSO  • 

0.0 

DFELC  « 

0.0 

DFELI  - 

0.0 

OFELO  * 

0.0 

DFEwC I » 

0.0 

DFEWCO* 

0.0 

DFEWI 1 « 

0.0 

DFEWI0= 

0.0 

DFEwOI = 

0.0 

DFEWOO" 

0.0 

DFXLE  * 

0.0 

DFXWI  - 

0.0 

DFXWO  « 

0.0 

DHB  ■ 

0.0 

DHCR  « 

0.0 

DHCT  - 

0.0 

DHTT  • 

0,0 

DISWl  « 

0.0 

DISWO  ■ 

0.0 

DLELC  « 

0.0 

DLEU  ■ 

0.0 

OLELO  = 

0.0 

OLESLC" 

0.0 

DLESL l 1 

0.0 

DLESLO- 

0.0 

DLEWC I = 

0.0 

DLEWC0= 

0 .0 

DLEWI I- 

0.0 

DLEWtO* 

0.0 

OLEWOI = 

0.0 

DLEWOO" 

0.0 

DCSWI  » 

0.0 

DOSWO  « 

0.0 

DVB  ■ 

0.0 

DVCR  ■ 

0.0 

over  - 

0.0 

DVTT  « 

0.0 

owe 

0.0 

OWCLE  « 

0.0 

DWCR  « 

0.0 

DWC  T ■ 

4.1 

DWTT  » 

0.416 

OWYENG* 

0.0 

F-1X1  WING  TIP 

.DEF-55H  SLAT  I NBD • EDGE  LOCATION  (FACTOR  OF  SEMI-SPAN) 

»DEF«SSH  SLAT  OUTBD.EOGE  LOCAT ION( FACTOR  OF  SEMI-SPAN) 

*DEF*5Sh  maxI«'jv  WING  LIFT  COEFFICIENT 
.DEFOSH  SURFACE  CHORO  TO  WING  CHORD  RATIO 

• D£F»  56  H RATIO  OF  INBOARD  END  OF  SPOILER  TO  WING  S£«t-SPAN 

•DEFsSSH  ratio  of  outboard  end  of  spoiler  TO  WING  semi-  pan 
.DEF«55h  Slat  EDGE  and  LEADING  EDGE  STATIONS  (FT) 

« DEF  * 55H  BODY  width  (FT! 

• DEF*5  5H  SURFACE  LENGTH  (IN) 

.OEFe56rt  SURFACE  CHORD  AT  INBOARD  END  (IN) 

•D£F«55h  SURFACE  CHORD  AT  OUTBOARD  END! IN) 

.D£F»55H  SURFACE  INBOARD  EDGE  LOCATION  (DECIMAL  OF  SEMI-SPAN) 
•DEF*55H  leading  EDGE  INBOARD  STATION  (FT) 

•OEF=55H  LEADING  EDGE  OUTBOARD  STATION  (FT) 

•DEF«55h  LENGTH  CF  LEADING  EDGE  OUTBOARD  CF  CENTER  SLAT  (FT) 
.DEF*S5H  LENGTH  OF  LEADING  EDGE  INBOARD  Or  CENTER  SLAT  (FT) 

.DEF «55H  LENGTH  Cc  FIXED  LEADING  EDGE  OUTBOARD  OF  SlAT  (FT) 

« CEF*55H  INBOARD  LEADING  EDGE  CHORD  OUTBOARD  OF  CENTER  SLAT-FT 
.DEFc55H  OUTBOARD  LEADING  EDGE  C-iGRD  OUTBOARD  OF  CENTER  SLAT-FT 
» DEF*55  h INBOARD  LEADING  EDGE  CHORD  INBOARD  CF  CENTER  SLAT-FT 

iDef«5Sh  outboard  leading  edge  c.-dad  j.<b^ard  of  center  slat-ft 

.DEF»55H  INBOARD  LEADING  EDor.  CHORD  OUTBOARD  CF  SLAT  (FT) 

«QEF ’55H  OUTSCARO  LEADING  EDGE  ChORO  OUTBOARD  OF  SLAT  (FT) 

.D-:F*55H  FIXED  LEADING  E JGE  LENGTH 

.DEF-55H  FIXED  LEADING  SOGE  wIDTh  INBD 

•DEF=S5h  FIXED  leading  EDGE  wIDTh  CUTED 

.DEF-55H  HORIZONTAL  TAIL  SPAN 

»DEF«55M  HO-  ZCNTAL  tail  ROOT  CHORD 

.OEF=53H  HORIZONTAL  TAIL  CHORD  AT  TIP  IFTI 

t DECr  5 5 w HORIZONTAL  TAIL  T-ICsNESS  AT  T I A (FT) 

.DEF-55M  INBOARD  SLAT  INBOARD  CHORD  (FEET) 

• DEF»5  5 H INBOARD  SLAT  OUTBOARD  CHORD  IFEET) 

• DEF  ^55H  LENGTH  CF  LEADING  EDGE  UNDER  CENTER  SLAT  i FT  ) 

.OSF-55H  LENGTH  OF  LEADING  EDGE  UNDER  INBOARD  SLAT  (FT) 

•DEF"S5H  length  of  leading  EDGE  UNDER  CUTbOARD  SLAT  i F T ) 

• DEF • 5 5 H CENTER  SLAT  LENGTH  (FEET) 

.DEF»55H  INBOARD  SLAT  LENGTH  (FEET) 

.DEF-55H  OUTBOARD  SLAT  LENGTH ( FEET > 

»OEF»55m  INBOARD  LEADING  EDGE  CHORD  UNDER  CENTER  SLAT  (FT) 
•DEF.55H  OUTBOARD  LEADING  EDGE  CHORD  UNDER  CENTER  SlAT  (FT) 

• DEF «5 5 H LEADING  EDGE  INBOARD  CHORD  UNDER  INBOARD  SLAT  (FT) 

• def.sbh  leading  e^ge  clt°oaru  chord  under  :n-*oakd  slat  t ft  i 
.OEF=55H  INBOARD  LEADING  EDGE  C«CRD  UNDER  OUTBOARD  SLA!  IFT) 
>OEF*55H  OUTBOARD  LEADING  EDGE  ChCRD  UNDER  OUTBOARD  SLAT  (FT) 
•DEF-55M  OUTBOARD  SLAT  INBOARD  CHORD  (FEET) 

»DEF»G5H  OUTBOARD  SLAT  OJT BOARD  ChORDIFEET) 

• DEF »55h  VERTICAL  TAIL  S = 4N 
*DEF»5SH  vertical  Tail  ROOT  CHORD 
•DEF«55h  VERTICAL  ChCRD  AT  T|P  |FT) 

• OEF.BiH  vertical  T h I C < N r S S A'  TIP  (FT) 

•DEF»55H  WING.HOSIZ  STAH.CR  v t R T STAB  SPAN  (FEET  I 

•DEFnSBH  wing  Chord  AT  SLAT  (FT!  actually  AT  MIDSPAN  OF  EXPOSED 
•OEF-5SH  WING  ChORO  AT  Th£  ROOT  (FT) 

• DEF « 5 5H  WING  CmORD  AT  TIP  (FT ) 

.DEF.55H  WING  ThICXNcSS  AT  ThE  TIP  (FT) 

•oefcssh  engine  station  (fti  two-engine  airplane 
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Table  4.  List  of  Input  (Continued). 


OWYIEN- 

0.0 

.OEF-55H 

DWYOEN- 

0.0 

•DEF-55H 

GHFS  - 

0.0 

•0EF-55H 

OHRS  - 

0.0 

i OEF-55  H 

GVFS  ■ 

0.0 

.DEF-55H 

GVRS  - 

0.0 

•0EF-55H 

GMAIL  - 

0.0 

,0EF=55H 

GWCFFL- 

0.0 

»0EF-=55h 

GWCFLP- 

0.0 

•DEF-55H 

GWCLE  - 

0.0 

.DEF-55H 

GWFS  - 

0.2A3A9 

.DEF.55H 

GVJRS  - 

0.69  166 

,D£F=55H 

GWTCR  - 

0.0 

• D E F • 5 5 H 

GWTCT  - 

0.0 

•0EF-55H 

NZ 

0.0 

.DEF-55H 

QENG  - 

0.0 

,0EF=55h 

QMAX  « 

0.0 

.DEF-55H 

SIGCR  - 

0.0 

,OEF -55H 

SW 

0.0 

.0EF-55H 

SWA  I L » 

0.0 

,DEF=55H 

SWFLAP- 

0.0 

•0EF-55H 

TIPWI  - 

18.  C 

•DEF-55H 

TSTl  - 

0.0 

.DEF-35H 

VS 

0.0 

•OEF-55H 

ZOLDGM- 

0.0 

.DEF-55H 

L 

0.0 

,0EF -55  H 

LSURF  « 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


SEND 


INBOARD  ENGINE  STATION  (FT)  FOUR-ENGINE  AIRPLANE 
OuTOOARD  ENGINE  STATION  ( FT  I FGUR-ENG I NE  AIRPLANE 
LOCATION  OF  HORIZONTAL  FRONT  SPAR  '(PERCENT  ROOT  CHORD) 
LOCATION  OF  HORIZONTAL  REAR  SPAR  (PERCENT  ROo T’  CnOnD  ) 
LOCATION  OF  VERTICA'.  FRONT  SPAR  (PERCENT  ROOT  CHORD) 
LOCATION  OF  VERTICA.  REAR  SPAR  (PERCENT  ROOT  CHORD) 
KATIO  OF  INBOARD  AIlEROW  STATION  TO  WING  SEMI-SPAN 
RATIO  OF  FOBEFLAP  CHORD  TO  FLAP  CHORD 
FLAP  CHORD  TO  wing  CHORD  RATIO 

ratio  of  slat  chord  to  wing  chord 

WING  wing  front  SPAR  LOCATION  (FACTOR  OF  WING  CHORD) 
LOCATION  OF  wing  REAR  SPAR  (.PERCENT  ROOT  CHORD) 

WING  I/C  AT  ROOT 

WING  T/C  AT  TIP 

DESIGN  LOAD  FACTOR 

NUMBER  OF  ENGINES 

MAXIMUM  DYNAMIC  PRESSURE  ( P S I ) 

DENSITY  FACTOR  AT  ALTITUDE 

wlNG.HORIZ.  STAB. .OR  VERT. STAB  AREA  (SQ.FT.) 


WING  FLAP  AREA  (SQ.FT. ) 
INPUT  TIP  WIDTH 
OESIGN  GROSS  WEIGHT  (LBS.) 
STALL  SPEED  KNOTS) 

FLAG  FOR  GECM  ROUTINES 
SURFACE  NU'-'3SR 
FOREFLAP-  1 THRU  20 
-21  THRU  2S 


7»DEF=55H 


LLI1 ) 

m 

0.0 

»DEF=55H 

LL  ( 2 ) 

m 

0.0 

i CEF  = 55  h 

LL  ( 3 ) 

n 

0.0 

.DEF-55 H 

ILIA) 

m 

0.0 

»DEF=55h 

LLI5I 

M 

0.0 

,0EF  = 55h 

LL  ( 6 ) 

* 

0.0 

»DEF=55h 

LL  ( 7 I 

* 

0.0 

• 0 E F ->  5 5 h 

LLiei 

n 

0.0 

*0EF  = 55H 

MATTY 

P*- 

2.0EF-53M 

FLAP 

AILERON  =26  ThRU  30 

slat  * a i thru  a& 

RUDDER  = A7  THRU  51 
ELEVAT0R=52  Thru  56 
SURFACE  CCDS 

1- AILERON 

2- RUOOER 

3- ELEVA  TOR 
A-FLAP 

5- SLAT 

6- FOREFLAP 

7- TIP 

8- FIXED  LEADING  EDGE 

9- SPOILERS  * 

10-FIXED  TRAILING  EDGE 

NO  SLAT S-EN TIRE  LENGTH  IS  FIXED  LEADING  EDGE 

INBOARD  SLA TS  ONLY 

CENTER  SI  - TS  ONLY 

outboard  slats  only 

I NBC ARO  AND  CENTER  SLATS  ONLY 

CENTER  A.Nu  CUTSCASD  SLATS  ONLY 

INBOARD  AND  OUTBOARD  SlATS  ONLY 

INBOARD, OUTBOARD  AND  CENTER  SLATS  ONLY 

MATERIAL  TYPE 

1 -7075  ALUMINUM 
-202A  ALUMINUM 
-DM5  17 9w  TITANIUM 
-17-7PH  STAINLESS  STEEL 
-5AE  1018  STEEL 

-BORON-EPOXY  COMPOS  I TE ( 0.005  IN) 

-BORON  ALUMINUM  COMPOSITE 
-GRAPHI TE-EPOXY  COMPOSITE 

-HEXCEL  ALUMINUM  hONE YCOM8  (1/8  HEX  0.001  FOIL) 


2 

3 

A 

5 

6 

7 

8 
9 
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iUi  k L •!  ♦ 'JuW/k'  '<£  « JJim  (iiSi 


1 


i 


II 


V: 


Table  5.  Aileron  Input. 


PSINPNL 

c 

C AX  WING  AILERONS 


MATTYP-2 4 
LSURF-1 . 
DW8-A7. 75 . 
DWCT-4 . 5 * 
GWFS-0.14, 
GWRS-0 . 5 3 5 t 
DWCR- 1 1.4* 
VS-100. 4 
QMAX-950. 4 
S I GCR “0 . ***«  8 6 ? 

MAX-2.5  . 
N2-7.5  4 
SW-380. 4 
TST 1-38300 . 4 
OCSL-72. 4 

ocsws - 19 . a 4 

DCSWO-13.54 

0CSY1-0.7A64 

SEND 


Table  7.  Elevator  Input. 


PS1NPNL 

c 

C C-5A  INNER  ELEVATOR  HORIZONTAL 


MATTYP-2  4 
LSURF • 3 4 
DWB-222 . 7 . 
DHCR-2 I.64 
OHCT-7.7  . 
GHR5-0 . 6 5 4 
VS-100.0, 
QMAX*  550. 4 
5 I GCR  -0  4 **“8  6 4 
CLMAX-4, 5 4 
NZ-3, 75  4 
SW-6200  • 0 4 
TST 1-769000 .0 1 
OCSL-193.35  4 

OCSWI  «8<*  4 0 4 

DCSWO-60  a 4 
SENO 


Table  6.  Rudder  Input. 


P8INPNL 

c 

C C-5A  RUDDER-LOWER 


LSURF -2  4 
MATTYP-2  4 
DWB-222. 7 4 
DVCR-30.93, 
DVCT*2<*  a 7 5 4 
GVRS-0.6A  4 
VS-100.04 
OMAX- 550 . 0 4 
5IGCR  = 0.a<.864 
CLMAX-h . 5 4 
NZ-3. 754 
SW-6200 . 0 4 
TST I- 769000 . 0^ 
DCSL-208 . 1 5 . 

DC SW 1-133. 63t 
DCSwO-133.634 
SEND 


Table  8.  Flap  Input. 


PSINPNL 

c 

C AX  WING  FLAPS 
C 

MATTYP-2  , 
LSURF-A , 

DWB-47, 75 4 
DWCT-A.5  4 
GWFS-0  a 1 A 4 
GWRS-0 . 5 3 5 . 
OWCR-ll.A, 
VS-100. 4 
0MAX-950 . 4 
S I GCRaO . a A86  a 
CLMAX- 2 . 5 4 
NZ- 7, 5 4 
SW-380. . 
TST1-383004  4 
DCSL-  182  4 4 
DCSWI  = 2 3 . , 
DCSWO-15. , 
GWCFLP-0.26  4 
SEND 


I 
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Table  9.  Slat  Input. 


PSINPNL 

C 

C AX  WING  SLATS 
C 

MATTYP-2. 
LSURF-5  > 
OWB-47.75. 

OWCT -4.5. 
GWFS-0.14, 
GWRS-0.53S. 
DWCR-11.4, 

VS-100.. 
QMAX-950. • 
SIGCR-0.4486. 
CLMAX*  2 » 5 » 
NZ-7.S, 

SW-380. * 
TST1-38300. . 
0(11-1.92. 
DI2I-5.0. 
0(31-6. 2» 
D(4)-17»0  . 
D(5)«17.5. 
0(61-23.0. 
OCSWI-1.3. 
0CSWO-0.8, 
0ISWI-1.6. 
DISWO-1.4. 
DLE5LC- 10.8 . 
DLESL 1=2.75 . 
OLESLO-5.83 . 
DOSWI -0.8 . 
OOSWO-O. 65 . 
CSCDWC-0  . 14** . 
*EnD 


Table  11.  Leading  Edge  Input 


PSINPNL 

C 

C F-lll  FIN  LEADING  EDGE 
C 

MATTYP-2 . 

OFXLE  -12.33  . 

OFXWI  -3.03  . 

OFXWO  -1.38  . 

L - 62. 

LSURF  - 8. 

SEND 


Table  10.  Fore  Flap  Input. 


PSINPNL 

c 

C AX  WING  FOREFLAPS 

c 

MATTYP-2. 

LSURF -6. 

6WCFFL-0 . 26 » 
GWCFLP-0.26 . 
OWB-47,75. 

DWCT-4.5 . 

GWFS-O. 14 . 

GWRS-0 . 535. 

DWCR*  1 1.4. 

VS-100.  . 

QMAX-9S0. » 
SIGCR-0.Ai.86. 

CLMAX*  2 . 6 . 

NZ-7.S. 

SW-380. . 

TST 1*38300.  . 
OCSL-182. . 

0CSWI-6. . 

DCSW0-4. » 

SEND 


Table  12.  Trailing  Edge  Input 


PSINPNL 

C 

C C-1A1  STAB  TE 
C 

MATTYP-2  . 

L-58. 

LSURF -10. 
GWTCR-O. 106 . 
GWTCT-0.106 . 
DHB-50 .36. 
GHRS-0 . 67  , 

OHCR* 1 a , 06  * 
DHCT-6.20. 
DBW-3.333  . 
0HB-50.35, 

.SEND 


Mm 


Table  13.  Spoilers  Input. 


P8INPNL 

c 

c f-111  SPOILERS 
C 


DWCT 

-3.824 

•DEF-55H 

GwRS 

■ 0.65 

.DEF-55H 

OWCR 

*12.506 

»D££:-55H 

DCSL 

-267.  A 

.DEF-55H 

DCSWI 

-21.7 

•OEc*55H 

DCSWO 

-11.8 

.DEF-55H 

LSURF 

a 

9 .DEF-55H 

mattyp 

“2  i 

CSCDwC-0.162 

• 

GWTCR 

-0.1135 

• 

GWTCT 

-0.085 

. 

GWCFLP 

-0.35 

. 

FLPTVP-0.0 

. 

DWU-60 

.533. 

ocsrl- 

SEND 

0.2561. 

LOCATION  OF  WING  REAR  SPAR  (PERCENT  ROOT  CHORD) 
WING  CHORD  AT  THE  ROOT  ( FJ ) 

SURFACE  LENGTH  (IN) 

SURFACE  CHORD  AT  INBOARD  END  (IN) 

SURFACE  Chord  AT  OUTBOARD  END  (IN) 

SURFACE  IDENTIFICATION  FLAG 


Table  14.  Sample  Output. 
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